SPECIFICATION 



TITLE OF THE INVENTION 
A COLOR TRANSFORMATION TABLE CREATING METHOD, A COLOR 
5 TRANSFORMATION TABLE CREATING APPARATUS, AND A 
COMPUTER READABLE RECORD MEDIUM IN WHICH A COLOR 
TRANSFORMATION TABLE CREATING PROGRAM IS RECORDED 

1) Field of the Invention 

10 The present invention relates to a method and 

apparatus for creating a color transformation table 
according to which color signals (for example, RGB 
signals) outputted from a color input device such 
as an image scanner or the like are correlated with 

15 color signals (for example, L*a*b* signals) in a color 
space which are different from a color space of the 
color input device, and a computer readable record 
medium in which a color transformation table creating 
program is recorded. 

20 2) Description of the Related Art 

Generally, a color space is a coordinate 
system expressing colors. There are various color 
spaces. A Color printer (color output device), an 
image scanner (color input device), or the like has 

25 a color space (device-dependent color space) 
depending on the device such as CMY (Cyan, Magenta 
and Yellow) or RGB (Red, Green and Blue) . A color 
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printer whose color signal values are from 0 to 255 
as the minimum value and the maximum value, for example, 
expresses output colors by color signal values of 
each of CMY in a range from 0 to 255. 
5 Other than the device-dependent color space, 

L*a*b* or XYZ is used, which is a color space 
(device-independent color space) not depending on 
a device. A color space not depending on a device 
can express absolute colors, thus is used as an 

10 intermediate color space between device-dependent 
color spaces when color transformation is performed 
between color spaces of different devices. Here, 
color transformation signifies that signal values 
expressed in a certain color space are transformed 

15 into color signal values expressed in another color 
space . 

For example, in order to match color 
appearance of a color original with color appearance 
of a printed matter of the color original obtained 
20 by reading it by a color scanner and outputting it 
from a color printer, color transformation is 
performed in the following procedures (1) and (2) : 

(1) Color- trans forming an RGB image obtained 
by the image scanner into an image expressed in, for 

25 example, the L*a*b* system. 

(2) Transforming the L*a*b* image obtained 
in (1) into CMY data for the color printer, and 
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outputting an image by the color printer on the basis 
of the CMY data. 

In color transformation according to the 
above procedures (1) and (2), color transformation 
5 tables (profiles) are used. In the color 

transformation table, relationships between 
different color spaces are registered. In color 
transformation, such the relationships are used. 
When a color registered in a color transformation 

10 table is transformed, the color transformation table 
is referred to read a relationship of the color 
registered therein, and color transformation is 
performed. When a color not registered in the color 
transformation table is transformed, an 

15 interpolating process is performed on the basis of 
data registered in the color transformation table 
to transform the desired color. 

A color transformation table (profile) 
showing input characteristics of a color input device 

20 such as an image scanner or the like is a collection 
of about 5000 relationships between color signal 
values (RGB values) showing coordinates in the RGB 
color system and chor omat icit y coordinate values 
(L*a*b* values) in the L*a*b* color space, as will 

25 be described later with reference to, for example, 
Table 1. 

When color transformation is performed using 
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such a color transformation table, it is important 
to save labor of creating a color transfo rma t ion table, 
and to create a color transformation table according 
to which accurate color transformation is possible. 
5 The present invention is aimed to create a 

color transformation table showing input 
characteristics of a color input device such as an 
image scanner. Hereinafter, description will be 
made of creation of a color transformation table for 

10 an image scanner. 

When a color transformation table is created, 
a standard color chart defined by ISO is read by an 
image scanner that is an object, then color signal 
values representing reflected light intensities of 

15 three primary color components are obtained from each 
color region in the color chart image. 

Next, chromaticity coordinate values 
representing an absolute color in each color region 
in the above standard color chart are measured by, 

20 for example, a colorimeter, the chromaticity 
coordinate values are correlated with color signal 
values corresponding to the relevant color region, 
whereby a relationship that becomes a base of a color 
transformation table is obtained. 

25 After that, an interpolating operation 

process is performed on the basis of the relationship 
that becomes the base, whereby a relationship between 
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color signal values corresponding to each of grid 
points uniformly distributed in the RGB space and 
chromaticity coordinate values of the absolute color 
in the L*a*b* color space that should correspond to 
5 the color signal values is determined. 

Namely, the measuring work as stated above 
is done on, for example, 288 colors expressed on the 
color patch chart, the interpolating operation 
process is performed on the basis of relationships 

10 of the 288 sets obtained through the measuring work, 
and relationships of 5000 sets are created. 

As a technique relating to the above 
interpolating operation process, there is a technique 
disclosed in, for example, Japanese Patent Laid-Open 

15 Publication No. 7-307872. The publication cites 
"Shepard' s Method" in which interpolation is 
performed with a weight inversely proportional to 
the fourth power of a distance, as a method for 
interpolation in a non-grid-shape color space. The 

20 method solves a problem of "scalloping" occuring in 
this method. 

Fig. 18 shows a simplified example of 
"scalloping". As shown in Fig. 18, when vectors A, 
B, C, D and E are measured, a new vector F between 

25 the vectors A and B is obtained by averaging all the 
vectors A, B, C, D and E weighted in the Shepard's 
Method. As a result, the interpolated vector F 
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differs from a vector indicated by a broken line 
estimated on the basis of the vectors A and B due 
to the contributions of the vectors C, D and E some 
distance away from the vector F. Such phenomenon is 
5 called "scalloping". Use of a color transformation 
table with narrow spaced grids is effective to 
appreciably avoid such "scalloping". 

In the technique disclosed in the above 
publication, the number of relationships configuring 
10 the color transformation table is increased by liner 
interpolation, the increased relationships are 
registered to create a color transformation table 
with a larger number of relationships (with shorter 
spacing) . The above technique can provide a color 
15 transformation table in which relationships larger 
in number than relationships actually measured, 
allowing labor saving in transformation table 
creation, preventing the effect of "scalloping", 
which results in accurate color transformation. 
20 On the other hand, there has been proposed 

a technique for estimating spectral characteristics 
of a color scanner on the basis of relationships 
obtained as results of measurement, and obtaining 
relationships configuring a color transformation 
25 table using the obtained spectral characteristics. 
As a technique for estimating input characteristics 
of a color input device such as an image scanner, 
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there is a technique disclosed in, for example, 
Japanese Patent Laid-Open Publication No. 6-46252. 
In the publication, there is proposed a technique 
for solving simultaneous equations derived on the 
5 basis of measured values in the least squares method 
to estimate the spectral characteristics. 

However, a work to measure a number of color 
samples and a work to read a number of color samples 
by a color input device that is an object are essential 
10 even when a color transformation table is created 
by interpolating basic relationships {Japanese 
Patent Laid-Open Publication No. 7-307872) , or when 
a color transformation table is created on the basis 
of spectral characteristics estimated from results 
15 of measurement (Japanese Patent Laid-Open 
Publication No. 6-46252). 

Particularly, the technique disclosed in 
Japanese Patent Laid-Open Publication No. 7-307872 
is effective to lighten a problem of "scalloping" 
20 by increasing the number of relationships by linear 
interpolation as stated above. However, the 
technique performs linear interpolation on original 
relationships, so that the number of relationships 
cannot be accurately increased if the number of 
25 relationships is originally small. In other words, 
it is essential to prepare a considerable number 
(several hundreds colors) of relationships that 
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become a base in order tocreateacolor transformation 
table attaining accurate color transformation. 
However, in order to prepare relationships that 
become the base of the color transformation table, 
5 it is essential to measure several hundreds of color 
signal values obtained by reading the original by 
a color input device and outputting it, and to measure 
several hundreds of color signal values in a color 
space which are different from a color space of the 
10 color input device, which is troublesome. 

Alternatively, there has been proposed a 
technique (disclosed in PCT Application No. 
PCT/ JP99/06447 ) , in which image color signal values 
(RGB values) obtained by reading a color chart by 
15 a color input device can be correlated with three 
weights on principal component vectors of spectral 
reflectance, respectively, by using estimated 
spectral characteristics of a color input device, 
whereby spectral reflectance or common color signal 
20 values (L*a*b* values, XYZ values or the like) of 
the color chart are calculated from image color signal 
values. According to this technique, the spectral 
characteristics can be estimated even if the number 
of color patches (the number of color regions in the 
25 color chart) is only several tens. As a result, a 
color transformation table can be createdwitha small 
number of color patches, leading to labor saving in 
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color transformation table creation. 

However, spectral reflectance of color 
patches ( color regions ) on the color chart has various 
shapes. For this, if only one kind of weights are 
5 used as weights on three principal component vectors 
of the spectral reflectance, that is, if the same 
color transformation formula is used for the whole 
color space to determine the spectral reflectance 
as a linear sum of three principal component vectors, 
10 color transformation accuracy in a part of the regions 
might largely degrade. As a result, it is impossible 
to accurately express the spectral reflectance, 
leading to degradation of accuracy of the color 
trans format ion table . 

15 

SUMMARY OF THE INVENTION 
In the light of the above problems, an object 
of the present invention is to provide a color 
transformation table creating method, a color 
20 transformation table creating apparatus and a 
computer readable record medium in which a color 
transformation table creating program is recorded, 
by which relationships of color signals are 
accurately obtained in a simplified way to create 
25 a color transformation table making accurate color 
transformation possible without much labor. 

The present invention therefore provides a method 
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for creating a color transformation table correlating 
a color signal outputted from a color input device 
in a color space of the color input device (hereinafter 
referred to as a transformation source color space) 
5 with a color signal in a color space (hereinafter 
referred to as a transformation target color space) 
which is different from the transformation source 
color space, comprising the steps of, a dividing step 
of beforehand defining a plurality of regions 

10 obtained by dividing the whole of the transformation 
target color space, and a color transformation table 
creating step of creating the color transformation 
table by using a plurality of color transformation 
formulas corresponding to the plural regions, 

15 respectively, defined at the dividing step. 

The method may further comprise a reading step 
of reading a plurality of color regions on a color 
chart by the color input device, and outputting a 
color signal in the transformation target color space 

20 corresponding to each of the respective color regions 
from the color input device, a colorimetric step of 
measuring the plurality of color regions by a 
colorimeter, and outputting spectral reflectance 
corresponding to each of the color regions from the 

25 colorimeter, a classifying step of classifying the 
spectral reflectance according to which region among 
the plural regions in the transformation target color 
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space a color signal in the transformation target 
color space corresponding to the spectral reflectance 
belongs to, and a spectral characteristics estimating 
step of estimating spectral characteristics of the 
5 color input device on the basis of the color signal 
outputted from the color input device at the reading 
step and the spectral reflectance outputted from the 
colorimeter at said colorimetric step, wherein, at 
the color transformation table creating step, the 
10 color transformation formula is created for each of 
the regions in the transformation target color space 
on the basis of the spectral reflectance classified 
at the classifying step and the spectral 
characteristics estimated at the spectral 
15 characteristics estimating step. 

The color transformation table creating step 
may comprise a relationship creating step of creating 
a relationship between a color signal in the 
transformation source color space and a color signal 
20 in the transformation target color space by using 
the color transformation formula according to each 
region, for each of the color transform formulas, 
and a creation processing step of creating the color 
transformation table on the basis of the relationship 
25 created for each of the color transformation formulas 
at the relationship creating step. 

The present invention further provides an 
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apparatus for creating a color transformation table 
comprising a color transformation table creation unit 
for creating the color transformation table by using 
a plurality of color transformation formulas 
5 corresponding to a plurality of regions, respectively, 
the regions being obtained by dividing the 
transformation target color space. 

The present still further provides a computer 
readable record medium in which a color 
10 transformation table creating program for making a 
computer realize a function of creating a color 
transformation table is recorded, the color 
transformation table creating program making the 
computer function as a color transformation table 
15 creation unit for creating the color transformation 
table by using a plurality of color transformation 
formulas corresponding to a plurality of regions, 
respectively, the regions being obtained by dividing 
the transformation target color space. 
20 According to the color transformation table 

creating method, the color transformation table 
creating apparatus and the computer readable record 
medium in which the color transformation table 
creating program is recorded, relationships all over 
25 the color space can be accurately calculated, by using 
color transformation formulas (color transformation 
formulas differing from divided region to divided 
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region) according to respective divided regions in 
the transformation target color space when a color 
transformation table for a color input device is 
created . 

5 Accordingly, it is possible to create a color 

transformation table attaining accurate color 
transformation with a small number of color patch 
data without much labor. By using a color 
transformation table created as above, it is possible 

10 to improve color transformation accuracy all over 
the color space. 

Principal component vectors to spectral 
reflectance calculated for each divided region have 
a high rate of contribution to the spectral 

15 reflectance existing in a corresponding divided 
region. Spectral reflectance obtain by measuring 
colors in plural kinds of color regions on the color 
chart is classified according to divided region, and 
a relationship is calculated by using a color 

20 transformation formulas corresponding to the divided 
region on the basis of the classified spectral 
reflectance and estimated spectral characteristics 
of the color input device. It is thereby possible 
to calculate a relationship more accurately than when 

25 principal component vectors to spectral reflectance 
calculated from the whole color space are used. 
According to this invention, principal component 
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vectors are calculated from spectral reflectance of 
a color region (color patch) included in each divided 
region. Since spectral reflectance of two colors at 
a short distance from each other in a color space 
5 is similar, it becomes easy to express spectral 
reflectance with principal component vectors, thus 
a relationship can be calculated more accurately, 
as a result . 

A plurality of relationships are created by 
10 using a plurality of color transformation formulas 
with respect to one color signal, and a color 
transformation table is created on the basis of the 
plural relationships. It is therefore possible to 
calculate a more accurate relationship, and create 
15 a color transformation table which attains more 
accurate color transformation. Using a color 
transformation table created as above, it is possible 
to further improve color transformation accuracy all 
over the color space. 
20 The present invention provides the following 

functions and effects: 

[1] As the transformation target color space, 
a uniform color space, in which a difference between 
colors when viewed by man and a distance between the 
25 colors in the color space are almost equal (CIELAB, 
CIELUV, CIECAM) , is used, whereby a distance in the 
color space of any color can be used as a weighting 
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factor . 

[2] Since a plurality of regions obtained by 
dividing the transformation target color space have 
regions overlapping on each other, gradation of color 
5 transformation result is assured. Namely, a result 
of interpolation operation on plural color 
transformation results calculated by using color 
transformation formulas for respective divided 
regions are used as a color transformation result 
10 existing in the overlapped region, therefore 
gradation of a color in each region is assured. 

[3] Use of plural regions obtained by dividing 
the whole transformation target color space according 
to hue angle makes accurate calculation of 
15 relationships possible. Namely, when principal 
component vectors are calculated on the basis of a 
spectral reflectance group obtained by classifying 
spectral reflectance of a color chart into hue angle 
regions, the principal component vectors have a high 
20 rate of contribution to the spectral reflectance 
group belonging to a corresponding hue angle region. 
Accordingly, such the principal component vectors 
can express a spectral reflectance group belonging 
to a corresponding hue angle region more accurately, 
25 thus a relationship can be accurately calculated, 
as a result. Such effect is remarkable when a color 
chart including colors having various hue angles is 
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used . 

[4] A hue angle of a primary color (for example, 
Cyan, Magenta, Yellow or the like) is the center of 
each region, whereby relationships can be calculated 

5 accurately. A shape of spectral reflectance of a 
color existing in the vicinity of the center of a 
primary color is similar to a shape of spectral 
reflectance of the primary color. For this, when 
principal component vectors are calculated on the 

10 basis of a spectral reflectance group existing in 
the vicinity of the center of a primary color, a rate 
of contribution of the principal component vectors 
to the spectral reflectance group in the vicinity 
of the center of the primary color is high, as a matter 

15 of course, thus a relationship can be calculated 

accurately . 

[5] Use ofplural regions obtained by dividing 
the whole transformation destination color space 
according to chroma makes accurate calculation of 

20 relationships possible. Namely, when principal 
component vectors are calculated on the basis of a 
spectral reflectance group obtained by classifying 
spectral reflectance of a color chart into chroma 
regions, a rate of contribution of the principal 

25 component vectors to the spectral reflectance group 
belonging to a corresponding chroma region is high, 
as a matter of course. Accordingly, such the 
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principal component vectors can express the spectral 
reflectance group belonging to its chroma region more 
accurately, thus a relationship can be accurately 
calculated, as a result. Such effect is remarkable 
5 when a color chart including colors having various 
kinds of chroma is used. 

[6] Use of plural regions obtained by dividing 
the whole transformation tar get color space according 
to lightness makes accurate calculation of 

10 relationships possible. Namely, when principal 
component vectors are calculated on the basis of a 
spectral reflectance group obtained by classifying 
spectral reflectance of a color chart into lightness 
regions, a rate of contribution of the principal 

15 component vectors to the spectral reflectance group 
belonging to a corresponding lightness region is high, 
as a matter of course. Accordingly, such the 
principal component vectors can express the spectral 
reflectance group belonging to its lightness region 

20 more accurately, thus a relationship can be 
accurately calculated, as a result. Such effect is 
remarkable when a color chart including colors having 
various kinds of lightness is used. 

[7] When a color transformation result, into 

25 which one color signal in the transformation target 
color space is transformed through a color 
transformation formula, belongs to a region 
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corresponding to the color transformation formula, 
the color transformation result is determined to be 
correct, and a color transformation table is created 
on the basis of the color transformation result 
5 determined to be correct. It is therefore possible 
to accurately calculate a relationship. Color 
transformation accuracy of a color transformation 
formula obtained from a corresponding divided region 
in a color space is low when a color signal outside 
10 region, to which the color transformation formula 
belongs, is color-transformed. For this, by 
discriminating a color transformation result having 
low accuracy from a color transformation result 
having high accuracy, it is possible to obtain a color 
15 transformation result having higher accuracy. 

[8] When there are a plurality of color 
transformation results with respect to one color 
signal determined to be correct, a color 
transformation result with respect to the color 
20 signal is calculated on the basis of values relating 
to distances between the plural color transformation 
results determined to be correct and boundaries 
(boundary surfaces) of regions to which the 
respective color transformation results belong, 
25 therefore gradation of the color transformation is 
assured. There is a case where the plural color 
transformation results calculated by using 
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respective color transformation formulas largely 
differ from one another and shift from the respective 
region boundary surfaces. In such case, gradation 
can be assured by performing a calculating process 
5 (interpolating process) on color transformation 
results as above. 

[9] When there are a plurality of color 
transformation results with respect to one color 
signal determined to be correct, one of the plural 
10 color transformation results determined to be correct 
is selected as a color transformation result with 
respect to the color signal on the basis of values 
relating to distances between the plural color 
transformation results determined to be correct and 
15 boundaries (boundary surfaces) of regions to which 
the color transformation results belong. It is 
thereby possible to shorten the processing time 
required for relationship creation. When the number 
of divided regions is large, a time required for the 
20 calculating process ( interpolating proces s ) on color 
transformation results becomes long. In order to 
shorten the processing time, it is determined that 
one at the largest distance from the boundary surface 
among correct color transformation results has the 
25 highest reliability among the plural correct color 
transformation results, and that color 

transformation result is determined as a final color 
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trans format ion result. 

[10] When there is no color transformation 
result determined to be correct with respect to one 
color signal, a color transformation result with 
respect to the color signal is calculated on the basis 
of reciprocals of values relating to distances 
between the plural color transformation results 
obtained with respect to the color signal and 
boundaries (boundary surfaces) of regions to which 
the respective color transformation results belong. 
Even if all the color transformation results are 
incorrect, it is possible to calculate a color 
transformation result having higher reliability. 
Namely, when all the color transformation results 
are incorrect, one at the smallest distance from the 
boundary surface among the incorrect results is 
determined to have higher reliability. 

[11] When there is no color transformation 
result determined to be correct with respect to one 
color signal, one of plural color transformation 
results obtained with respect to the color signal 
is selected as a color transformation result with 
respect to the color signal on the basis of reciprocals 
of values relating to distances between the plural 
color transformation results obtained with respect 
to the color signal and boundaries (boundary 
surfaces) of regions to which the respective color 
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transformation results belong. Even if all the color 
transformation results are incorrect, it is possible 
to obtain a color transformation result having higher 
reliability within a short time. 
5 [12] The color transformation table is used 

to correlate a color signal in the transformation 
source color space with spectral reflectance 
according to a color transformation result as a color 
signal in the transformation target color space. It 
10 is thereby possible to express a color transformation 
result by spectral reflectance, and substitute the 
color input device for a simplified colorimeter. 

[13] When there are a plurality of color 
transformation results determined to be correct with 
15 respect to one color signal, one kind of spectral 
reflectance of the color signal is calculated on the 
basis of values relating to distances between the 
plural color transformation results determined to 
be correct and boundaries (boundary surfaces) of 
20 regions to which the respective color transformation 
results belong, whereby gradation of the spectral 
reflectance is assured. 

[14] When there are a plurality of color 
transformation results determined to be correct with 
25 respect to one color signal, one kind of spectral 
reflectance according to the plural color 
transformation results determined to be correct is 
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selected as spectral reflectance with respect to the 
color signal on the basis of values relating to 
distances between the plural color transformation 
results determined to be correct and boundaries 
5 (boundary surfaces) of regions to which the 
respective color transformation results belong. It 
is thereby possible to shorten a processing time 
required for relationship creation. 

[15] When there is no color transformation 
10 result determined to be correct with respect to one 
color signal, spectral reflectance with respect to 
the color signal is calculated on the basis of 
reciprocals of values relating to distances between 
plural color transformation results obtained with 
15 respect to the color signal and boundaries (boundary 
surfaces) of regions to which the respective color 
transformation results belong. Even if all the color 
transformation results are incorrect, it is possible 
to calculate spectral reflectance having higher 
20 reliability. 

[16] When there is no color transformation 
result determined to be correct with respect to one 
color signal, one kind of spectral reflectance 
according to plural color transformation results is 
25 selected as spectral reflectance with respect to the 
color signal on the basis of reciprocals of values 
relating to distances between the plural color 
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transformation results obtained with respect to the 
color signal and boundaries (boundary surfaces) of 
regions to which the respective color transformation 
results belong. Even if all the color transformation 
5 results are incorrect, it is possible to obtain 
spectral reflectance having higher reliability 
within a short time. 

[17] As values relating to distances between 
plural color transformation results and boundaries 
10 (boundary surfaces) of regions to which the 
respective color transformation results belong, 
valuesof AE94 obtained by transforming the distances 
are used, whereby a color transformation table 
attaining more accurate color transformation when 
15 viewed by man is created. Since AE94 has a color 
difference closer to man's view than a distance (A 
Eab) in the L*a*b* space, a color transformation 
result using AE94 is more accurate when viewed by 
man . 

20 

BRIEF DESCRIPTION OF THE DRAWINGS 
Fig. 1 is a block diagram showing a functional 
structure of a color transformation table creating 
apparatus according to an embodiment of this 

25 invention; 

Fig. 2 is a block diagram showing a hardware 
structure of a computer system (color transformation 
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table creating apparatus) for accomplishing a color 
transformation table creating method according to 
the embodiment of this invention; 

Fig. 3 is a diagram showing relationships 
5 between the RGB space and the L*a*b* space; 

Figs. 4A through 4C are diagrams for 
illustrating an example of region division according 
to hue angle according to the embodiment; 

Fig. 5 is a flowchart for illustrating a flow 
10 of the color transformation table creating method 
according to the embodiment of this invention; 

Fig. 6 is a flowchart for illustrating 
procedures for estimating spectral characteristics 
according to the embodiment; 
15 Fig. 7 is a flowchart for illustrating 

procedures for creating a color transformation 
formula and procedures for creating a relationship 
according to the embodiment; 

Fig. 8 is a flowchart for illustrating 
20 procedures for determination and procedures for an 
interpolating process according to the embodiment; 

Figs. 9 through 12 are diagrams for 
illustrating a practical example of the interpolating 
process according to the embodiment; 
25 Fig. 13 is a diagram for illustrating an 

example of region division according to chroma 
according to a first modification of the embodiment; 
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Figs. 14 through 16 are diagrams for 
illustrating a practical example of the interpolating 
process according to the first modification of the 
embodiment ; 

5 Fig. 17 is a diagram for illustrating an 

example of region division according to lightness 
according to a second modification of the embodiment; 
and 

Fig. 18 is a diagram showing a simplified 
10 example of scalloping. 

DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Hereinafter, description will be made of an 
embodiment of the present invention. 

15 [1] Description of a structure of an embodiment 
Fig. 1 is a block diagram showing a functional 
structure of a color transformation table creating 
apparatus according to an embodiment of this 
invention. Fig. 2 is a block diagram showing a 

20 hardware structure of a computer system (color 
transformation table creating apparatus) for 
accomplishing a color transformation table creating 
method according to the embodiment of this invention. 
Fig. 3 is a diagram showing relationships between 

25 the RGB space and the L*a*b* space. Figs. 4A through 
4C are diagrams (a*b* plane) for illustrating an 
example of region division according to hue angle 
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according to the embodiment. 

The computer system (personal computer: PC) 

1 according to this embodiment shown in Figs. 1 and 

2 functions as an apparatus for creating a color 
5 transformation table representing relationships 

between a color space (transformation source color 
space: the RGB space in this embodiment) of an image 
scanner 60 as a color input device and a color space 
(transformation target color space: the L*a*b* space 

10 in this embodiment ) independent of a device . Bymeans 
of the color transformation table, it is possible 
to correlate RGB signals outputted from the image 
scanner with signals in the L*a*b* space. 

The color transformation table for a color 

15 input device generally stores relationships in terms 
of about 5000 (to 30000) colors. In this embodiment, 
description will be made by way of creation of a 
grid-like color transformation table in the RGB space 
(RGB space having 17x17x17 grid points) having values 

20 in 17 levels in each of RGB colors, as shown in Table 
1 below. Here, ranges of the L*a*b* values are -100 
^ L* ^ 100 , -100 ^ a* ^ 100 , and -100 ^ b* ^ 100, for 
example . 
[Table 1] 
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L*a*b*values 


(0,0,0) 

(0,0,15) 

• 
• 
• 

(255,255,255) 


( L *(0,o,o) a *(o,o,o) *<o,o,o>) 

l 1 * *(0,0,15> a *(0,0A5) b *(0,0,15)) 
• 

• 

• 

( L *(255,255,255) a *<255,255,255) * (255,2 5 5, 255) ) 



The grid-like color transformation table in 
the RGB space represents relationships between signal 
values of colors (colors in grid points) distributed 
5 in grid arrangement in the RGB space and color signal 
values of the colors in the L*a*b* space, that is, 
relationships between RGB values obtained by reading 
a color original by the image scanner and outputting 
it, and L*a*b* colorimet ric values of the color 

10 original. The color transformation table stores 
"grid-like RGB value data sets" in which values 
regularly increase as RGB values as shown in, for 
example, Table 1. 

RGB values are distributed in grid 

15 arrangement in the RGB space, but the RGB space does 
not have a linear relationship with the L*a*b* space. 
For this, points in the L*a*b* space obtained by 
transforming RGB values into L*a*b* values do not 
distribute in grid arrangement , distribution of which 

20 is distorted. Fig. 3 conceptually shows 
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relationships between the RGB space and the L*a*b* 
space in two dimensions. 

This embodiment is described by way of example 
where the transformation target color space is a 
5 uniform color space (CIELAB, CIELUV, CIECAM or the 
like) , in which a difference between colors viewed 
by man is almost equal to a distance between the colors 
in the color space, particularly, the L*a*b* space. 

The L*a*b* space is a color space defined by 

10 CIE (Commission Internationale de 1 ' Eclair age ) , in 
which a distance in the space represents a difference 
in color appearance that man senses. Namely, it is 
possible to use a distance in the L*a*b* space as 
a weighting factor for color transformation on any 

15 color in the L*a*b* space. To the contrary, in the 
CIE 1931 standard colorimetric system (xy 
chromaticity diagram) or the Munsell color system, 
a distance and a difference in color appearance do 
not correspond, thus a distance in the color space 

20 cannot be used as a weighting factor for color 
transformation. However, there are other coordinate 
systems such as the L*u*v* space or the CIECAM space 
other than the L*a*b* space, in which a distance and 
a difference in color appearance are equal, like the 

25 L*a*b* space. It is possible to use them. 

As shown in Fig. 2, the computer system 1 is 
provided with a CPU 10 and a bus line 50 connected 
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to the CPU 10. The bus line 50 is connected to a memory 
(ROM or RAM) 20, along with a display (for example, 
CRT, LCD, PDP or the like) 30, a keyboard 40 and a 
mouse 41 via an input/output interface 51. 
5 A state of display on the display 30 is 

controlled by the CPU 10 to display various 
information for use in creation of a color 
transformation table. The keyboard 40 and the mouse 
41 are operated by the operator (user) who refers 

10 to a screen of the display 30 to input various 
directions and various information to the CPU 10 
(computer system 1) . 

The memory 20 stores an application program 
(color transformation table creating program) 21 for 

15 realizing a color transformation table creating 
apparatus, which comprises a spectral 

characteristics estimation unit 11, a spectral 
reflectance classification unit 12, a color 
transformation table creation unit 14 and an RGB space 

20 grid point signal value output unit 15 (shown in Fig. 
1) . 

The CPU 10 reads the above application program 
21 from the memory 20 to execute it, thereby 
accomplishing functions as the spectral 
25 characteristics estimation unit 11, the spectral 
reflectance classification unit 12, the color 
transformation table creation unit 14 and the RGB 
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space grid point signal value output unit 15 (details 
of which will be described later) , that is, a function 
of the color transformation table creating apparatus 
of this invention. As a result, a color 

5 transformation table creating method of this 
invention is executed. 

The above application program 21 is recorded 
in a computer readable recordmedium such as a flexible 
disk, a CD-ROM or the like, and provided. The 
10 computer system 1 reads the program 21 from the record 
medium, transfers and stores it in an internal storage 
unit or an external storage unit to use the same. 
Alternatively, the program 21 may be recorded in a 
storage unit (record medium) such as a magnetic disk, 
15 an optical disk, a magneto-optic disk or the like, 
and provided to the computer system 1 from the storage 
unit through a communication path. 

When the functions as the spectral 
characteristics estimation unit 11, the spectral 
20 reflectance classification unit 12, the color 
transformation table creation unit 14 and the RGB 
space grid point signal value output unit 15 are 
attained, the program 21 stored in the internal 
storage unit (memory 20 in this embodiment) is 
25 executed by a microprocessor (CPU 10 in this 
embodiment) of the computer. At this time, the 
computer system 1 may directly read the program 
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recorded in a record medium to execute it. 

In this embodiment, a computer is a concept 
including hardware and an operation system, 
signifying hardware operating under control of the 
5 operation system. When the operation system is 
unnecessary and the application program solely 
operates the hardware, the hardware itself 
corresponds to the computer. The hardware comprises 
at least a microprocessor such as a CPU or the like, 
10 and a means for reading the computer program recorded 
in a record medium. 

The application program 21 includes program 
codes for realizing the functions as the spectral 
characteristics estimation unit 11, the spectral 
15 reflectance classification unit 12, the color 
transformation table creation unit 14 and the RGB 
space grid point signal value output unit 15 by such 
the computer (computer system 1). A part of the 
functions may be realized by an operation system, 
20 not by the application program 21. 

As a record medium in this embodiment, any 
one of various computer readable record media such 
as an IC card, a ROM cartridge, a magnetic tape, a 
punched card, an internal storage unit (memory such 
25 as an RAM or an ROM) of the computer, an external 
storage unit, a printed matter on which codes such 
as bar codes are printed, and the like may be employed, 
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other than the above flexible disk, CD-ROM , DVD, 
magnetic disk, optical disk and magneto-optic disk. 

In the computer system 1 shown in Figs. 1 and 
2, an image scanner (color input device) 60 and a 
5 colorimeter 70 are connected to the bus line 50 via 
an input/output interface 52. 

The computer system 1 of this embodiment 
creates a color transformation table for the image 
scanner 60 on the basis of results of reading of a 
10 color chart 100 (refer to Fig. 1) by the image scanner 
60, and results of colorimetry of the color chart 
100 by the colorimeter 70. 

The input/output interface 52 functions as 
an input unit for inputting RGB signals corresponding 
15 to color regions obtained by reading several tens 
of color regions on the color chart 100 by the image 
scanner 60, and results of colorimetry (spectral 
reflectance and L*a*b* values to be described later) 
of the color regions by the colorimeter 70. The RGB 
20 signals and the results of colorimetry inputted 
through the input /output interface 52 are temporarily 
stored in the memory 20. 

The colorimeter 70 is used to measure each 
color region on the color chart when a color 
25 transformation table for a color input device is 
created, which outputs spectral reflectance of 
380-730 nm at resolution of 10 nm. The spectral 
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reflectance outputted at this time is normalized with 
the maximum value of 1.0. 

The colorimeter 70 of this embodiment has a 
function of transforming spectral reflectance into 
5 L*a*b* values and outputting the values, along with 
a function of measuring and outputting the spectral 
reflectance. The L*a*b* values are required when the 
CPU 10 (spectral reflectance classification unit 12) 
classifies the spectral reflectance, as will be 

10 described later. 

According to this embodiment, L*a*b* values 
corresponding to spectral reflectance are obtained 
by the colorimeter 70. Alternatively, the CPU 10 
(spectral reflectance classification unit 12) may 

15 transform spectral reflectance from the colorimeter 
70, and calculate it. 

The color chart 100 has at least several tens 
of color patches (color regions) . As the color chart 
100, an input color target ( AN SI/IT8 . 7/2) by ANSI 

20 (American National Standards Institutes) standard 
may be used. 

In such case, the input color target has a 
color patch region including 264 color patches, and 
a gray patch region (gray scale) including gray 

25 patches in different 24 gradation levels. In this 
embodiment, suitable several tens of color regions 
among these color regions are designated, and the 
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designated color regions are read by the image scanner 
60, or measured by the colorimeter 70. 

The memory 20 also stores a color 
transformation table 22 created by the CPU 10, along 
5 with the above application program (color 
transformation table creating program) 21. 

The memory 20 functions as a working memory 
when the color transformation table 22 is created, 
as will be described later with reference to Fig. 
10 1. The memory 20 stores RGB signals from the image 
scanner 60, spectral reflectance (L*a*b* values ) from 
the colorimeter 70, spectral characteristics, 
classification information, relationships 1 to 3, 
and the like . 

15 Next, description will be made in more detail 

of various functions (functions as the spectral 
characteristics estimation unit 11, the spectral 
reflectance classification unit 12, the color 
transformation table creation unit 14 and the RGB 

20 space grid point signal value output unit 15) 
accomplished by the CPU 10. 

The spectral characteristics estimation unit 
11 reads out RGB signals and spectral reflectance 
of the same color chart 100 (the same color patch) 

25 obtained by the image scanner 60 and the colorimeter 
70, respectively, from the memory 20, and estimates 
spectral characteristics of the image scanner 60 on 
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the basis of the RGB values and the spectral 
reflectance. Estimation of the spectral 

characteristics is done in a method that will be 
described later with reference to Expressions (1) 
5 to (10) . 

The spectral reflectance classification unit 
12 determines which region in the L*a*b* space L*a*b* 
values corresponding to each spectral reflectance 
obtained by the colorimeter 70 belong to, and 

10 classifies the spectral reflectance of several tens 
of color patches obtained by the colorimeter 70 into 
groups according to belonging region. Results of the 
classification are outputted as classification 
information to the memory 20, and stored therein. 

15 According to this embodiment, three regions 

1, 2 and 3 obtained by dividing the whole L*a*b* space 
are beforehand defined and set (dividing step). 
These three regions (hereinafter referred to divided 
regions, occasionally) 1 to 3 are obtained by dividing 

20 the L*a*b* space by boundaries (practically boundary 
surfaces in parallel to the L* axis) Bl, B2 and B3 
according to hue angle, as shown in Figs. 4A through 
4C, for example. 

As shown in Fig. 4A, the region 1 is a region 

25 in a range at an angle of 120° centered at a hue angle 
of 210° of the primary color C (Cyan). As shown in 
Fig. 4B, the region 2 is a region in a range at an 
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angle of 120° centered at a hue angle of 330° of the 
primary color M (Magenta) . As shown in Fig. 4C, the 
region 3 is a region in a range at an angle of 120° 
centered at a hue angle of 90° of a primary color Y 
5 (Yellow) . 

A region having a width of d (constant) as 
shown in Figs. 4A through 4C is added to and formed 
in each of the regions 1 to 3 in order to overlap 
the regions 1 to 3 on one another. By forming such 

10 a region, overlapped regions as shown in Figs. 9 
through 12, that is, a region formed by overlapping 
the regions 1 to 3 on one another, a region having 
a width of 2d formed by overlapping the regions 1 
and 2 on each other, a region having a width of 2d 

15 formed by overlapping the regions 2 and 3 on each 
other, and a region having a width of 2d formed by 
overlapping the regions 3 and 1 on each other, are 
formed . 

The color transformation table creation unit 
20 14 forms the color transformation table 22 using color 
transformation formulas corresponding to the 
respective above-described three regions 1 to 3, 
which comprises a relationship creation unit 141, 
a determination unit 142 and a creating process unit 
25 143 . The RGB space grid point signal value output 
unit 15 generates signal values (grid point signal 
values) of colors distributed in grid arrangement 
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in the RGB space as described above with reference 
to Table 1, and outputs the signal values to the color 
transformation table creation unit 14 (relationship 
creation unit 141) . 
5 The relationship creation unit 141 creates 

color transform formulas for respective regions in 
the L*a*b* space on the basis of spectral reflectance 
groups classified by the spectral reflectance 
classification unit 12 and spectral characteristics 

10 of the image scanner 60 estimated by the spectral 
characteristics estimation unit 12, and creates 
relationships between RGB values (color signals in 
the transformation source color space) from the RGB 
space grid point signal value output unit 15 and L*a*b* 

15 values (color signals in the transformation target 
color space) , by using the three color transformation 
formulas . 

Incidentally, the color transformation 
formulas are created in a method to be described later 

20 with reference to Expressions (11) to (16), 
accompanying creation of the relationship by the 
relationship creation unit 141. A color 

transformation formula used to transform RGB values 
into L*a*b* values is practically created by 

25 combining Expressions (14) to (16) with a formula 
used to transform XYZ values into L*a*b* values (not 
described since it is a known general formula) . 
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At this time, the relationship creation unit 
141 calculates three L*a*b* values of each of all 
grid points in the RGB space using three color 
transformation formulas to create three kinds of 
5 relationships, outputs the three kinds of 
relationships as relationships 1 to 3 to the memory 
20, and stores them in the memory 20. 

The determination unit 142 determines whether 
a result of color transformation (L*a*b* values) of 

10 RGB signals ( signal values of one grid point ) obtained 
by using a color transformation formula when the 
relationship creation unit 141 creates the 
relationships 1 to 3 belongs to a region corresponding 
to the color transformation formula. When the result 

15 belongs to the region, the determination unit 142 
determines that the result of color transformation 
is correct . 

In more detail, the determination unit 142 
determines whether the relationship 1 [L*a*b*l 

20 (Ll,al,bl) ] created by using a color transformation 
formula corresponding to the region 1 belongs to the 
region 1 shown in Fig. 4A. When the relationship 1 
belongs to the region 1, the determination unit 142 
determines that the relationship 1 is correct. 

25 Similarly, the determination unit 142 determines 
whether the relationship 2 [L*a*b*2 (L2,a2,b2)] 
created by using a color transformation formula 
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corresponding to the region 2 belongs to the region 
2 shown in Fig. 4B. When the relationship 2 belongs 
to the region 2 , the determination unit 142 determines 
that the relationship 2 is correct. The 
5 determination unit 142 also determines whether the 
relationships [L*a*b*3 (L3,a3,b3)] created by using 
a color transformation formula corresponding to the 
region 3 belongs to the region 3 shown in Fig. 4C. 
When the relationship 3 belongs to the region 3, the 

10 determination unit 142 determines that the 
relationship 3 is correct. 

The creating process unit 143 obtains a set 
of relationships on the basis of three relationships 
1 to 3 created using the three color transformation 

15 formulas by the relationship creation unit 141, and 
creates the color transformation table 22 in the 
memory 20. At this time, the creating process unit 
143 creates the color transformation table 22 on the 
basis of results of color transformation (L*a*b* 

20 values) determined to be correct by the determination 
unit 142. 

In more detail, when there are a plurality 
of results of color transformation determined to be 
correct by the determination unit 142 with respect 
25 to one grid point (one color signal) in the RGB space 
(when two or three results of color transformation 
are determined to be correct, in this embodiment), 
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the creating process unit 143 calculates a result 
of color transformation (L*a*b* values) of one grid 
point on the basis of distances between the plural 
results of color transformation (L*a*b* values) 
5 determined to be correct and boundaries of regions 
to which the respective results of color 
transformation belong, and outputs the result, as 
will be described later with reference to Figs. 9 
and 10, and Expressions (17) and (18). 

10 When there is one result of color 

transformation determined to be correct by the 
determination unit 142 with respect to one grid point 
in the RGB space, the creating process unit 143 outputs 
the result (L*a*b* values) of color transformation 

15 determined to be correct as it is as a result (L*a*b* 
values) of color transformation with respect to the 
grid point, as will be described later with reference 
to Fig. 11 and Expression (19). 

When there is no result of color 

20 transformation determined to be corrected by the 
determination unit with respect to one grid point 
in the RGB space, the creating process unit 143 
calculates a result of color transformation with 
respect to the grid point on the basis of reciprocals 

25 of distances between three results of color 
transformation obtained with respect to the grid 
point by the relationship creation unit 141 and 
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boundaries of regions to which the respective results 
of color transformation belong, and output it, as 
will be described later with reference to Fig. 12 
and Expression (20) . 
5 When there are a plurality of results of color 

transformation determined to be correct by the 
determination unit 142 with respect to one grid point 
in the RGB space, the creating process unit 143 may 
function to select a result of color transformation 

10 having the highest reliability, and output it. 

In such case, the determination unit 143 
selects one of a plurality of results of color 
transformation determined to be correct as a result 
of color transformation with respect to one grid point 

15 on the basis of distances between the plural results 
of color transformation determined to be correct and 
boundaries of regions to which the respective results 
of color transformation belong, and outputs it. At 
this time, the creating process unit 143 selects, 

20 for example, a result of color transformation having 
the largest distance {a result of color 
transformation at the largest distance from the 
boundary) . 

Similarly, when there is no result of color 
25 transformation determined to be correct by the 
determination unit 142 with respect to one grid point 
in the RGB space, the creating process unit 143 may 
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function to select a result of color transformation 
having the highest reliability, and output it. 

In such case, the creating process unit 143 
selects one of three results of color transformation 
5 as a result of color transformation with respect to 
the grid point on the basis of reciprocals of distances 
between the three results of color transformation 
obtained with respect to the grid point by the 
relationship creation unit 141 and boundaries of 

10 regions to which the respective results of color 
transformation belong, and output s it . At this time, 
the creating process unit 143 selects a result of 
color transformation having the largest reciprocal 
of the distance (a result of color transformation 

15 at the shortest distance from the boundary) , for 
example . 

In the example shown in Fig. 1, the 
relationships 1 to 3 created by the relationship 
creation unit 141 are temporarily stored in the memory 
20 20. Alternatively, the relationships 1 to 3 may be 
directly inputted to the determination unit 142 and 
the creating process unit 143 from the relationship 
creation unit 141. 
[2] Description of operations according to the 
25 embodiment 

[2-1] Description of procedures for creating a color 
t ran s format ion table 
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Next, description will be made of a flow of 
a color transformation table creating method 
(operations of the color transformation table 
creating apparatus 1 in the above structure) 
5 according to the embodiment of this invention, with 
reference to a flowchart (steps Sll to S25) shown 
in Fig . 5 - 

According to this embodiment , before creation 
of the color transformation table 22 is initiated, 

10 three regions 1, 2 and 3 obtained by dividing the 
whole L*a*b* space are beforehand defined and set 
(dividing step; not shown). 

When the color transformation table creating 
apparatus 1 of this embodiment creates the color 

15 transformation table 22 for the image scanner 60, 
as shown in Fig. 5, several tens of color patches 
on the color chart 100 are read by the image scanner 
60 (step Sll; reading step), RGB signals 
corresponding to each of the color patches are 

20 outputted to the color transformation table creating 
apparatus 1 from the image scanner 60, and stored 
in the memory 20 through the input/output interface 
52 (step S12; reading step) . 

The same several tens of color patches as those 

25 read by the image scanner 60 are measured by the 
colorimeter 70 (step S13; colorimetric step), 
spectral reflectance corresponding to each of the 



43 



color patches is outputted from the colorimeter 70, 
and stored in the memory 20 through the input/output 
interface 52 (step S14; colorimetric step) . At this 
time, the colorimeter 70 of this embodiment 
5 transforms the measured spectral reflectance into 
L*a*b* values corresponding to each of the color 
patches, and the L*a*b* values are stored in the memory 
20 through the input/output interface 52 (step S15) . 

The spectral characteristics estimation unit 

10 11 estimates spectral characteristics of the image 
scanner 60 on the basis of the RGB signals and the 
spectral reflectance of each of the color patches 
stored in the memory 20 according to procedures to 
be described later with reference to Expressions (1) 

15 to (10) (step S16; spectral characteristics 
estimating step) . The spectral characteristics are 
outputted to the memory 20 and stored therein (step 
S17) . 

The spectral reflectance classification unit 
20 12 determines which region 1, 2 or 3 shown in Fig. 
4A, 4B or 4C the L*a*b* values of each of the color 
patches belong to. The spectral reflectance of each 
of the color patches is classified into a group 1, 
2, ... or K (K=3 in this embodiment) according to the 
25 region (step S18; classifying step) . A result of the 
classification is stored as classification 
information in the memory 20 (step S19). 
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After that, the color transformation table 
creation unit 14 creates the color transformation 
table 22 while creating three color transformation 
formulas corresponding to the respective regions 1 
5 to 3 in the L*a*b* space on the basis of the spectral 
reflectance groups (groups 1 to 3) classified by the 
spectral reflectance classification unit 11 and the 
spectral characteristics of the image scanner 60 
estimated by the spectral characteristics estimation 
10 unit 12 (steps S21 to S25; color transformation table 
creating step ) . 

At this time, the RGB space grid point signal 
value output unit 15 generates RGB values (grid point 
signal values) of each of colors distributed in grid 
15 arrangement in the RGB space, and outputs the RGB 
values to the color transformation table creation 
unit 14 (relationship creation unit 141) (step S20) . 

The relationship creation unit 141 creates 
relationships between the RGB values (grid point 
20 signal values) at step S21 and L*a*b* values by each 
of the color transformation formulas, by using the 
three color transformation formulas corresponding 
to the respective region 1 to 3 (stepS21; relationship 
creating step) . 
25 At this time, as color transformation 

formulas used to transform the RGB values into the 
L*a*b* values, Expressions (14) to (16) to be 
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described later and a formula used to transform XYZ 
values into L*a*b* values (known general formula) 
are employed. At step S21, L*a*b* values of each of 
all grid points in the RGB space are calculated by 
5 using the three color transformation formulas, and 
three kinds of relationships to 3 are created and 
outputted (step S22). 

Next, the determination unit 142 determines 
whether each set of L*a*b* values (a result of color 
10 transformation) included in the relationship 1, 2 
or 3 created at step S21 belong to the region 
corresponding to the color transformation formula, 
as stated above. When the L*a*b* values belong to 
the region, the determination unit 142 determines 
15 that the result of color transformation is correct, 
and outputs a result of the determination to the 
creating process unit 143 (step S23; determining 
step ) . 

The creating process unit 143 determines a 
20 relationship on the basis of the relationships 1 to 
3 created using the color transformation formulas 
at step S21, thereby creating the color 
transformation table 22 (steps S24 and S25) . At this 
time, basically on the basis of L*a*b* values 
25 determined to be correct at step 2 4 among L*a*b* values 
included in the relationships 1 to 3 with respect 
to a grid point in the RGB space, final L*a*b* values 
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of the grid point are determined (step S24; creating 
process step) . The L*a*b* values determined as this 
are outputted as one establishing the relationship 
with respect to the grid point to the memory 2 0, whereby 
5 the color transformation table 22 is created in the 
memory 20 (step S25) . 

When two or three results of color 
transformation are determined to be correct, the 
creating process unit 143 calculates or selects a 

10 result of color transformation (L*a*b* values) of 
the grid point on the basis of distances between the 
two or three results of color transformation (L*a*b* 
values) determined to be correct and boundaries of 
regions to which the respective results of color 

15 transformation belong, and determines it. 

When a result of color transformation 
determined to be correct is only one, the creating 
process unit 143 outputs the result of color 
transformation (L*a*b* values) determined to be 

20 correct as it is as a result of color transformation 
(L*a*b* values) of the grid point, as stated above. 

To the contrary, when no result of color 
transformation determined to be correct exists, the 
creating process unit 143 calculates or selects a 

25 result of color transformation of the grid point on 
the basis of reciprocals 1 of distances between three 
results of color transformation (L*a*b* values) 
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obtained with respect to the grid point by the 
relationship creation unit 141 and boundaries of 
regions to which the respective results of color 
transformation belong, and determines it, as stated 
5 above. 

[2-2] Description of procedures for spectral 
characteristics estimation 

Next, detailed description will be made of 
a process at step S16 in Fig. 5, that is, procedures 
10 for spectral characteristics estimation (operation 
of the spectral characteristics estimation unit 11) 
according to this embodiment, with reference to a 
flowchart (steps S31 to S35) shown in Fig. 6. 

The spectral characteristics estimation unit 
15 11 (step S16) estimates the spectral characteristics 
on the basis of several tens of color signals and 
several tens of kinds of spectral reflectance in the 
following procedures. 

[2-2-1] Creation of simultaneous linear formulas 
20 (step S31) 

Here, the spectral characteristics are 
treated as discrete data corresponding to wavelengths 
(380, 390, 730 nm) for every 10 nm in a visible 

radiation wavelength range 380-730 nm. 
25 Generally, a value obtained by multiplying 

the spectral characteristics by spectral reflectance 
of a color patch is RGB values. 
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In order to estimate spectral characteristics 
of the image scanner 60, at step S31, simultaneous 
linear formulas for each of R, G and B comprising 
36 or more equations with spectral characteristics 
being unknown are created as shown in Expression (1) 
below, using RGB values obtained by reading color 
patches by the image scanner 60 and spectral 
reflectance of the color patches obtained by the 
colorimeter 70. 
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In the above Expression (1), R±, Gi and Bi (i 
= 1, ... N> are values (RGB values) obtained by reading 
N color patches by the image scanner 60, RE F { i , j ) (i 

15 = 1, N; j = 380, 390, 730) is N sets of spectral 
reflectance data of the N color patches obtained by 
the colorimeter 70, and REF { i, 38 o}/ REF ( i 7 390 }/ 
REF(i,730 ) are obtained spectral reflectance (36 
discrete data) of the i-th color patch. LS Rj (j = 

20 380, 390, 730) is spectral characteristics for 

R (36 discrete data), LS Gj (j = 380, 390 , 730) is 

spectral characteristics for G (36 discrete data), 
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and L S b j (j = 380, 390 , 730) is spectral 

characteristics for B (36 discrete data). 

When the number of color patches N is less than 
36, 36 simultaneous linear equations cannot be 
created. In such case, at step S31, simultaneous 
linear equations comprising 36 linear equations are 
created as shown by Expression (2) below. 
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[2-2-2] Singular value decomposition of 
simultaneous linear equations (steps S32 and S33) 
Here, a part of the above Expression (1) or 
(2) [that is, Expression (3) below] is called a 
spectral reflectance data set. 
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At step S32, singular value decomposition is 
performed on the spectral reflectance data set. 
Afterthat, atstepS33, noise included in the spectral 
reflectance data is removed on the basis of a result 



50 



of the singular value decomposition. 

Noise removal is performed in the following 
procedures. First, singular value decomposition is 
performed on the spectral reflectance data set, so 
that the spectral reflectance data set shown in the 
above Expression (3) i s decomposed into three met rics 
as shown by Expression (4) below. 
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In the above Expression (4), REF ( i,j)(i = 1, 
10 ... N, j = 380, 390, 730) is spectral reflectance 

data of N color patches obtained by the colorimeter 
70, as stated above, and W m (m = 1, 36) is weights 
by which principal component vectors are to be 
multiplied. P(m,j) (m = 1/ 36; j = 380, 390, 

15 730) is 36 sets of principal component vectors for 
the spectral reflectance. Each set of principal 
component vectors is 36 discrete data. U im (i = 1, 
N; m = 1 , 36) is coefficients. 

Weights Wi-W 36 calculated in the singular 
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value decomposition represent weights (degrees of 
importance) on the principal component vectors P{i,j) 
-P{36,j) ( j = 380, 390 , 730) . In proportion to each 
of the weights Wi-W 36/ a degree of importance of the 
5 principal component vectors corresponding thereto 
increases . 

Considering significant digits of the 
spectral reflectance data set, a weight at a distance 
of equal to or more than the significant digits from 

10 the largest weight can be assumed to be meaningless 
principal component vectors. At step S33, such 
principal component vectors are eliminated to remove 
noise included in the spectral reflectance data. 

Namely, at step S33, a threshold value 

15 (normally called a condition number) is determined 
by using Expression (5) below, and ones equal to or 
less than the threshold value among the weights Wi-W 36 
are replaced with "0" as shown in Expression (6) below 
to remove the noise. 

20 Threshold value (condition number) = [a 

maximum value among the weights Wi~W 36 ] / [significant 
digits of the data set (colorimeter)] ... (5) 

Incidentally, Expression (6) below is used 
when there are M weights that exceed the above 

25 threshold value, wherein the M weights are Wi-W M , and 
weights W M+1 -W 36 on the (M+l)th to 36th principal 
component vectors are replaced with w 0" . Where, M<36, 
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as a matter of course. 
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[2-2-3] Calculation of spectral characteristics 
(steps S34 and S35) 
5 Here, the matrices shown in the above 

Expression (6) are replaced with matrices as shown 
by Expression (7) below, for the sake of simplicity. 
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Where, U is a N*36 matrix having the above 
10 coefficients U im (i = 1, N; m = 1, 36) as elements, 
and V is a 36x36 matrix having the principal component 
vectors P (Itl/j) (m = 1, 36; j = 380 , 390, 730) 

as elements. The matrices U and V in the above 
Expression (7) have characteristics as shown in 
15 Expressions (8-1) and (8-2) below. Incidentally, in 
the Expressions (8-1) and (8-2), I is a unit matrix, 
U T is a transposed matrix of the matrix U, and V T is 
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transposed matrix of the matrix V, 
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The above Expression (1) or (2) is expanded 
as shown by Expression (9) below by using Expressions 
(8-1) and (8-2). By using Expression (10) below 
obtained from the expansion (9), spectral 
characteristics of the image scanner 60 is calculated 
(step S34), and the calculated spectral 
characteristics are outputted to the memory 20 and 
stored therein (step S35). 
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[2-3] Description of procedures for spectral 
reflectance classification 

Next, description will be made of a process 
at step S18 in Fig. 5, that is, procedures for spectral 
reflectance classification (operation of the 
spectral reflectance classification unit 12) 
according to this embodiment , with reference to Figs. 
4A through 4C. 

According to this embodiment, hue angles of 
the primary colors C, M and Y of the color chart 100 
are set to 210°, 330° and 90°, and the CIELAB space 
(L*a*b* space) is beforehand divided into three 
regions 1 to 3 with the above hue angles 210° , 330° 
and 90° being the centers, respectively, as shown 
in Figs. 4A through 4C. In order to form overlapped 
regions, a region having a width of d is added to 
and formed in each of the regions 1 to 3, as stated 
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above. The width d is set to, for example, 5.0. 

The spectral reflectance classification unit 
12 (step S18) determines which region 1, 2 or 3 shown 
in Fig. 4A, 4B or 4C L*a*b* values of each color patch 
5 belong to. Whereby, spectral reflectance of each of 
color patches is classified into any one of three 
groups 1 to 3 according to region. 

Namely, spectral reflectance whose L*a*b* 
values are determined to belong to the region 1 is 

10 classified into the group 1. Spectral reflectance 
whose L*a*b* values are determined to belong to the 
region 2 is classified into the group 2. Spectral 
reflectance whose L*a*b* values are determined to 
belong to the region 3 is classified into the group 

15 3 . 

According to this embodiment, L*a*b* values 
used for determination when the regions are divided 
are obtained by the colorimeter 70. Alternatively, 
the spectral reflectance classification unit 12 may 
20 transform spectral reflectance from the colorimeter 
70 into L*a*b* values corresponding to the spectral 
reflectance of each color patch. 

[2-4] Descriptions of procedures for color 
transform formula creation and procedures for 
25 relationship creation 

Next, detailed description will be made of 
processes at steps S20 to S22 in Fig. 5, that is, 
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procedures for color transform formula creation and 
procedures for relationship creation (operations of 
the relationship creation unit 141) according to this 
embodiment, with reference to a flowchart (steps S41 
5 to S47) shown in Fig. 7. 

The relationship creation unit 141 calculates 
L*a*b* values corresponding to grid-like RGB values. 

Namely, the relationship 1 all over the color 
space is calculated on the basis of grid-like RGB 
10 values, the spectral characteristics of the image 
scanner 60, and the principal component vectors 
calculated from spectral reflectance of the group 

1. Similarly, the relationship 2 all over the color 
space is calculated on the basis of grid-like RGB 

15 values, the spectral characteristics of the image 
scanner 60, and the principal component vectors 
calculated from spectral reflectance of the group 

2. The relationship 3 all over the color space is 
calculated on the basis of grid-like RGB values, the 

20 spectral characteristics of the image scanner 60, 
and principal component vectors calculated from 
spectral reflectance of the group 3. Whereby, three 
sets of L*a*b* values (L*a*b*l, L*a*b*2 and L*a*b*3) 
are calculated as results of color transformation, 

25 by using three kinds of color transformation formulas 
for one set of grid-like RGB values. 

The grid-like RGB values are outputted from 



57 



the RGB space grid point signal value output unit 
15. The spectral characteristics of the image 
scanner 60 are LS Rj (j = 380, 390, 730), LS Gj (j = 

380, 390, 730) and LS Bj (j - 380, 390, 730) 

5 estimated and calculated by using the above 
Expression (10) by the spectral characteristics 
estimation unit 11. The principal component vectors 
of each of the groups 1 to 3 are calculated as P(m,j) (m 
= 1, M; j = 380, 390, 730) using the above 

10 Expressions (4) to (6) by performing singular value 
decomposition on the spectral reflectance of each 
of the groups 1 to 3. 

Hereinafter, description will be made of 
procedures for calculating a relationship of a 

15 certain group (spectral reflectance group) with 
reference to a flowchart (steps S41 to S47) shown 
in Fig . 7 . 

[2-4-1] Calculation of weights by which principal 
component vectors are to be multiplied (steps S41 

20 and S42) . 

Spectral reflectance is reproduced as a linear 
sum of M sets of principal component vectors 
calculated by performing singular value decomposing 
method on the spectral reflectance groups (groups 

25 1, 2 and 3 in this embodiment) , as shown Expression 
(11) below. Incidentally, Expression (11) is 
obtained by extracting elements in the i-th row in 
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the above Expressions (4) and (6), 



(REF (i r 380) REF (i ,390) REF (i . 7 30)) 



= ( Un • Wi ) • ( P (i, 380) P(l, 390) 



P ( 1 , 730) ) 



+ ( U i2 " W 2 ) * ( P (2, 380 ) P{2 , 3 90) 



P (2, 730) ) 



5 



+ + (UiM # W M ) • (P<M, 3 8 0 ) P(M r 390) 



P (M, 7 3 0) ) 



(ID 



Where, REF (i/j) ( j - 380, 390, 730) is the i-th 
spectral reflectance classified into each of the 
groups 1 to 3. Incidentally, i is 1, N, and N is 
10 the number of spectral reflectance REF (i , j} classified 
into each of the groups 1 to 3. P< m/ j)(m == 1, M; 
j = 380 , 390, 730) is M sets of principal component 
vectors. U im -W m (m = 1, M) is weighting factors 

corresponding to the i-th spectral reflectance, by 
15 each of which the m-th principal component vectors 
(P ( m, 380), P<m,39o>, ... , P<m,73o>) are to be multiplied. 

As shown in the above Expression (11), it is 
possible to express spectral reflectance that is 36 
dimensional data by a weighting factor in about 
20 several dimension (3-7 dimension; M = 3-7) , by using 
the principal component vectors. 

RGB values (RGB) obtained by reading a certain 
color patch by the image scanner 60 and outputting 
it from the image scanner 60 are calculated on the 
25 basis of spectral characteristics LS Rj (j = 380, 390, 
730) , L S g j (j - 380, 390, 730), LS Bj (j = 380, 390, 
730) and spectral reflectance RE Fj (j = 380, 390, 
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730) obtained by measuring a color original of 
the color patch, as shown by Expression (12) below. 
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The spectral reflectance REFj ( j = 380, 390, 
730) is expressed as a linear combination of the M 
sets of principal component vectors as shown in the 
above Expression (11), Expression (12) thus can be 
rewritten into Expression (13) below. 
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A 

If inverse matrices of matrices indicated by 
A in the above Expression (13) can be calculated, 
weighting factors U im -W m (m = 1, M) are calculated 
from the RGB values as shown by Expression (14) below. 
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Hereinafter, this embodiment will be described 
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by way of example where M=3 (that is, the number of 
principal component vectors is three) . In this case, 
the matrices indicated by A become a 3><3 matrix. By 
substituting the RGB values (RGB) into the above 
5 Expression (14), weighting factors Ui-W x , U 2 -W 2 and 
U 3 -W 3 corresponding to the RGB values are uniquely 
calculated . 

According to this embodiment, singular value 
decomposition is performed on spectral reflectance 

10 of each of the groups 1 to 3 f whereby three sets of 
principal component vectors of each of the groups 
1 to 3 are calculated by using the above Expressions 
(4) to (6) (Step S41) • 

On the basis of the calculated principal 

15 component vectors P (mj) (m = 1, 2, 3; j = 380, 390 , 
730) and spectral characteristics LS Rj (j = 380, 
390, 730), LS Gj (j = 380, 390, 730) and LS Bj (j 

= 380, 390, 730), three kinds of Expressions (14) 
are created. By substituting RGB values from the RGB 

20 space grid point signal value output unit 15 in the 
Expression (14), weighting factors Ui-Wi, U 2 -W 2 and 
U 3 -W 3 by which the three sets of principal component 
vectors are to be multiplied are calculated (step 
S42) . 

25 [2-4-2] Calculation of spectral reflectance 

Next, the three sets of principal component 
vectors P ( i,j)(j = 380, 390 , 730), P (2 ,j>(j = 380 , 
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390, 730) and P ( 3,j>(j = 380 , 390 , 730) are 

multiplied by the weighting factors Ui-Wx, U 2 -W 2 and 
U 3 -W 3 obtained for each RGB values, respectively, by 
using the above Expression (14), and a linear sum 
5 thereof is determined. In consequence, spectral 
reflectance RERj (j = 380 , 390, 730) corresponding 
to the RGB values is calculated as shown by Expression 
(15) below (step S43) . 

(REF 38 o REF390 REF730) 

10 = (Ui ■ Wi) •( P (1 , 380 ) P(l,390> P<1,730)) 

+ ( U 2 * W 2 ) * ( P {2 , 380 ) P<2 , 390) P(2, 730)) 

+ (U 3 *W 3 ) • (P(3 , 380 ) P(3, 390) P(3, 730)) »■ (15) 

At this time, the spectral reflectance is 
corrected in such a way that when a value of the 

15 spectral reflectance RERj (j - 380, 390, 730) 
calculated using the above Expression (15) is less 
than 1.0 (that is, when the value is negative), the 
value is rewritten to 0.0, and when the value exceeds 
1.0, the value is rewritten to 1 . 0 . (stepS44). Using 

20 the corrected spectral reflectance, XYZ values are 
calculated using Expression (16) below (step S45) . 
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Where Lj • Xj ( j = 380, 730) is the product of color 
matching function X and light source, Lj • Yj ( j = 380, 
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730) is the product of color matching function 

Y and light source, and Lj-Zj (j = 380, 730) is 

the product of color matching function Z and light 

source . 

5 After the XYZ values are calculated from the 

RGB values as above, the XYZ values calculated by 
using the above Expression (16) are transformed into 
L*a*b* values by using the known transformation 
formula (step S46) . A relationship between the RGB 
10 values and the L*a*b* values in grid arrangement is 
thus created. 

The above relationship is created for each 
of the groups 1 to 3 that are spectral reflectance 
groups. As a result, three relationships 1 to 3 are 
15 created for each set of the RGB signal values (grid 
point signal values) as shown in, for example, Table 
2 below, and outputted to and stored in the memory 
20 (step S47 ) . 

According to this embodiment, the color 
20 transformation formula for transforming RGB values 
into L*a*b* values is a combination of the above 
Expressions (14) to (16) and the formula for 
transforming XYZ values into L*a*b* values, as stated 
above. The color transformation table creation unit 
25 14 (relationship creation unit 141) creates the 
relationships 1 to 3 between RGB values and L*a*b* 
values, while creating the above color transformation 
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formula for each of classified spectral reflectance 
groups, that is, for each divided region. 

[2-5] Description of procedures for determination 
and procedures for a creating process 
5 Next, processes at step S23 and S24 in Fig. 

5, that is, procedures for determination and 
procedures for a creating process (operations of the 
determination unit 142 and the creation process unit 
143) according to this embodiment will be described 

10 in detail with reference to Figs. 9 through 12 and 
a flowchart ( steps S51 to S58 ) shown in Fig. 8. Figs. 
9 through 12 are diagrams for illustrating a practical 
example of an interpolating process by the creation 
process unit 143. 

15 The creation process unit 143 performs an 

interpolating operation process (interpolating 
process) on three sets of L*a*b* values created with 
respect to one set of RGB values by the relationship 
creationunit 141 to calculate one set of L*a*b* values, 

20 for example. 

Prior to the calculating process 
(interpolating process) , the determination unit 142 
determines whether a result (L*a*b* values) of color 
transformation of an RGB signal obtained by the 

25 relationship creation unit 141 belongs to a 
corresponding region, according to this embodiment. 
When the result belongs to the corresponding region, 
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it is determined that the result of color 
transformation is correct. 

An example of three relationships 1 to 3 
created by the relationship creation unit 141 and 
5 an example of results of determination by the 
determination unit 142 are shown in Table 2 below. 
InTable 2, "O" is enteredina column of determination 
result when the determination unit 142 determines 
that the result of color transformation (L*a*b* 
10 values) is correct, while W 2C" is entered in a column 
of determination result when the determination unit 
142 determines that the result of color 
transformation (L*a*b* values) is incorrect. 
[Table 2] 
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As shown in Figs, 9 through 12, for example, 
a result of color transformation of the relationship 
1 is L*a*b*l (Ll,al,bl), and a distance between 
5 L*a*b*l (LI , al , bl ) and a boundary surface Bl of the 
region 1 is Dl . Similarly, a result of color 
transformation of the relationship 2 is L*a*b*2 
(L2,a2,b2), and a distance between L*a*b*2 (L2,a2,b2) 
and a boundary surface B2 of the region 2 is D2 . A 
10 result of color transformation of the relationship 
3 is L*a*b*3 (L3,a3,b3), and a distance between 
L*a*b*3 (L3,a3,b3) and a boundary surface B3 of the 
region 3 is D3 . 
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At this time, the determination unit 142 
calculates the distance Dl between L*a*b*l (LI, al, bl) 
and the boundary surface Bl of the region 1* The 
determination unit 142 adds 1 to AX iNumber_Of __True" 
5 that is a counter value when it is determined that 
L*a*b*l (Ll,al,bl) is inside the region 1, that is, 
it is determined that L*a*b*l is correct (step S51) . 

Similarly, the determination unit 142 
calculates the distance D2betweenL*a*b*2 (L2,a2,b2) 
10 and the boundary surface B2 of the region 2. When 
it is determined that L*a*b*2 (L2,a2,b2) is inside 
the region 2, that is, it is determined that L*a*b*2 
(L2,a2,b2) is correct, the determination unit 142 
adds 1 to "iNumber__Of__True" as a counter value (step 
15 S52) . The determination unit 142 further calculates 
the distance D3 between L*a*b*3 (L3,a3,b3) and the 
boundary surface B3 of the region 3. When it is 
determined that L*a*b*3 (L3,a3,b3) is inside the 
region 3, that is, it is determined that L*a*b*3 
20 (L3a3b3) is correct, the determination unit 142 adds 
1 to "iNumber__Of_True" as a counter value (step S53) . 

The creation process unit 143 calculates and 
determines one color transformation result L*a*b*4 
(L4,a4,b4) according to a value of the result of 
25 determination by the determination unit 142 (a value 
of "iNumber_Of_True) , that is, the number of " O " in 
Table 2) using Interpolation Expressions (17) to (20) 
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below (steps S54 to S58). 

When three color transformation results 
L*a*b*l (Ll,al,bl), L*a*b*2 (L2,a2,b2) and L*a*b*3 
(L3,a3,b3) belong to the regions 1 to 3, respectively 
5 as shown in, for example, Fig. 9, and 
"iNumber_Of_Ture" = 3 [refer to color transformation 
results of RGB values (0,0,0) inTable2, for example ] , 
the creation process unit 143 calculates and 
determines one L*a*b*4 (L4,a4,b4) on the basis of 
10 the distances Dl to D3 by using interpolation formulas 
as shown by Expressions (17) below (step S54). 
L4=(L1*D1+L2*D2+L3*D3)/(D1+D2+D3) 
a4=(al*Dl + a2*D2-ha3*D3)/ (D1+D2+D3) ... (17) 
b4= (bl*Dl+b2*D2+b3*D3) / (D1+D2+D3) 
15 Whentwocolor trans format ion results L*a*b* 1 

(Ll,al,bl) and L*a*b*2 (L2,a2,b2) belong to the 
regions 1 and 2 , respectively, as shown in, for example, 
Fig. 10, and "iNumber_Of_Ture" = 2 [refer to color 
transformation result of RGB values (0,0,15) in Table 
20 2, for example], oneL*a*b*4 (L4,a4,b4) is calculated 
and determined on the basis of the distances Dl and 
D2 by using interpolation formulas as shown by 
Expressions (18) below. 

L4=(L1*D1+L2*D2)/(D1+D2) 
25 a4= (al*Dl + a2*D2) / (D1 + D2) ... (18) 

b4=(bl*Dl+b2*D2)/(Dl+D2) 

When only one color transformation result 
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L*a*b*l (Ll,al,bl) belongs to the region 1 as shown 
in, for example, Fig. 11, and w iNumber_Of_True ,; = 
1 [refer to color transformation results of RGB values 
(0,0,31) inTable2, forexample], L*a*b*l (Ll,al,bl) 
5 is determined as it is as L*a*b*4 (L4,a4,b4) as shown 
by Expressions (19) below (step S56). 
L4=L1 

a4 = al ... (19) 
b4 = bl 

10 When all three color transformation results 
L*a*b*l (Ll,al,bl), L*a*b*2 (L2,a2,b2) and L*a*b*3 
(L3,a3,b3) do not belong to the regions 1 to 3, 
respectively, as shown in, for example, Fig. 12, and 
M i Number Of True " = 0 [refer to color transformation 

15 results of RGB values (255, 255, 255) in Table 2, for 
example], one L*a*b*4 (L4,a4,b4) is calculated and 
determined on the basis reciprocals of the distances 
Dl to D3 by using interpolation formula shown by 
Expressions (20) below (step S57). 

20 L4= (LI* ( 1/D1) +L2* (1/D2) +L3* (1/D3))/((1/D1 

) + (1/D2) + (1/D3) ) 

a4=(al*{l/Dl)+a2*(l/D2)+a3*(l/D3))/((l/Dl 

) + (1/D2) + (1/D3) 

b4=(bl*(l/Dl) +b2* (l/D2)+b3*(l/D3))/((l/Dl 
25 ) + (1/D2) + (1/D3) ) ... (20) 

L^a^b* values (L*a*b*4 (L4,a4,b4)) after the 
interpolating process obtained in any one of the 
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processes at steps S54 to S57 are outputted as a final 
result of color transformation corresponding to one 
set of RGB values to the memory 20, and a relationship 
of the RGB values is determined. 
5 The above processes are carried out on all 

grid-like RGB values (grid point signal values), 
whereby the color transformation table 22 for the 
image scanner 60 is created in the memory 20. 

When there are a plurality of color 

10 transformation results that are determined to be 
correct, the above creation process unit 143 
calculates and determines one L*a*b*4 (L4,a4,b4) on 
the basis of the distances Dl to D3 as shown in the 
Expressions (17) or (18). Alternatively, the 

15 creation process unit 143 may determine that a color 
transformation result having the largest distance 
(color transformation result at the largest distance 
from the boundary) has the highest reliability, and 
select and determine it as L*a*b*4 (L4,a4,b4). 

20 In the example shown in Fig. 9, L*a*b*l 

(Ll,al,bl) is selected and determined as a final color 
transformation result L*a*b*4 (L4,a4,b4) in the 
selecting process since Dl>D3>D2. In the example 
shown in Fig. 10, L*a*b*l (Ll,al,bl) is selected and 

25 determined as a final color transformation result 
since D1>D2 . 

Similarly, when there is no color 
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transformation result that is determined to be 
correct, the creating process unit 143 calculates 
and determines one L*a*b*4 (L4,a4,b4) on the basis 
of reciprocals of the distances Dl to D3 as shown 
5 in Expressions (20). Alternatively, the creating 
process unit 143 may determine that a color 
transformation result having the largest reciprocal 
of the distance (that is, color transformation result 
having the shortest distance and closest to the 

10 boundary) has the highest reliability, and select 
and determine it. 

In the example shown in Fig. 12, L*a*b*3 
(L3,a3,b3) is selected and determined as a final color 
transformation result L*a*b*4 (L4,a4,b4) in the 

15 selecting process since D3<DKD2. 

[3] Description of functions and effects of the 
embodiment 

According to this embodiment, it is possible 
to accurately calculate relationships all over the 

20 color space by using color transformation formulas 
according to respective divided regions of the L*a*b* 
space (that is, color transformation formulas 
differing from divided region to divided region) when 
the color transformation table 22 for the image 

25 scanner 60 is created. 

The relationship creation unit 141 of this 
embodiment calculates principal component vectors 
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for each divided region on the basis of a spectral 
reflectance group of color patches belonging to each 
divided region beforehand set. Two colors at a short 
distance from each other in the color space have 
5 similar spectral reflectance. For this, principal 
component vectors calculated on the basis of a 
spectral reflectance group belonging to the same 
region has a high rate of contribution to spectral 
reflectance existing in each divided region. Use of 

10 such principal component vectors are useful to 
express spectral reflectance, which in turn allows 
a relationship to be calculated highly accurately, 
as a result. Namely, Use of not principal component 
vectors calculated from the whole space but principal 

15 component vectors calculated fromeach divided region 
allows calculation of more accurate relationships. 

Accordingly, it is possible to create the 
color transformation table 22 attaining accurate 
color transformation with a small number of color 

20 patches (color chart data: image color signals, 
spectral re flectance ) without much labor . Use of the 
color transformation table 22 created as this can 
improve accuracy of color transformation all over 
the color space. 

25 The color transformation table creation unit 

14 creates three relationships 1 to 3 with respect 
to one RGB signal (grid-like signal values) by using 
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three color transformation formulas, and creates the 
color transformation table 22 on the basis of these 
relationships 1 to 3. It is therefore possible to 
calculate a relationship having higher accuracy in 
5 consideration of the relationships. As a result, it 
is possible to create the color transformation table 
22 attaining accurate color transformation, and 
improve the color transformation accuracy all over 
the color space by using the color transformation 
10 table 22 created as this. 

This embodiment provides the following 
effects . 

According to this embodiment, by using the 
L*a*b* space (CIELAB) as the transformation target 

15 color space, which is a uniform color space in which 
a difference between colors viewed by man and a 
distance between the colors in the color space are 
approximately equal, it is possible to use a distance 
in the color space of any color in the color space 

20 as a weighting factor when color transformation is 
pe r formed . 

The three divided regions 1 to 3 obtained by 
dividing the L*a*b* space have regions overlapping 
on each other, whereby gradation of a color 
25 transformation result is assured. Namely, results 
of the interpolation operation on a plurality of color 
transformation results calculated through the color 
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transformation equations for respective divided 
regions are used as a color transformation result 
L*a*b*4 (L4,a4,b4) existing in the over lapped region, 
gradation of colors in the divided regions 1 to 3 
5 is assured. 

Use of the regions 1 to 3 obtained by dividing 
the L*a*b* space according to hue angle allows 
accurate calculation of relationships . Namely, when 
principal component vectors are calculated on the 

10 basis of a spectral reflectance group obtained by 
classifying spectral reflectance of a color chart 
into hue angle regions, the principal component 
vectors have, of course, a high rate of contribution 
to the spectral reflectance group belong to the 

15 cor responding hue angle region 1 , 2, or 3. Therefore, 
such principal component vectors can express more 
accurately a spectral reflectance group belonging 
to each of the hue angle groups 1 to 3, so that 
relationships can be calculated accurately as a 

20 result. Sucheffectis remarkably when a color patch 
including colors having various hue angles. 

Particularly, the L*a*b* space is divided 
into three regions 1 to 3 with hue angles of the primary 
colors (Cyan, Magenta and Yellow) being the centers 

25 of the respective regions. A shape of spectral 
reflectance existing in the vicinity of the center 
of a primary color is similar to the shape of spectral 
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reflectance of the primary color. Accordingly, when 
principal component vectors are calculated on the 
basis of a spectral reflectance group existing in 
the vicinity of the center of a primary color, the 
5 principal component vectors have, of course, a high 
rate of contribution to the spectral reflectance 
group in the vicinity of the primary color, so that 
relationships can be accurately calculated. 

When a color transformation result (L*a*b* 
10 values) obtained by transforming one set of RGB values 
(grid point signal values) through a color 
transformation formula belongs to a region 
corresponding to the color transformation formula, 
it is determined that the color transformation result 
15 is correct. The color transformation table 22 is 
created on the basis of the color transformation 
result determined to be correct, so that 
relationships can be more accurately calculated. 
Accuracy of color transformation of a color 
20 transformation formula obtained from a certain 
divided region of the color space is low when a color 
signal outside its own region is transformed, as a 
matter of course. Accordingly, by discriminating a 
color transformation result having low accuracy from 
25 a color transformation result having high accuracy, 
it is possible to obtain a color transformation result 
having higher accuracy. 
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Where there are a plurality of color 
trans-formation results (L*a*b* values) determined 
to be correct with respect to one set of RGB signals 
(grid point signal values), a color transformation 
5 result L*a*b*4 (L4 / a4 / b4) with respect to the RGB 
signals is calculated on the basis of distances Dl 
to D3 between the plurality (two or three) of color 
transformation results determined to be correct and 
boundaries of regions to which the respective color 

10 transformation results belong, so that gradation of 
the color transformation is assured. There is a case 
where plural color transformation results calculated 
by using respective color transformation formulas 
largely differ and they are shifted from respective 

15 region boundary surfaces. In such case, a 

calculating process ( interpolating process) as above 
is performed on such the color transformation results, 
so that gradation is assured. 

When there are a plurality of color 

20 transformation results (L*a*b* values) determined 
to be correct with respect to one RGB signal (grid 
point signal values) , one of the plural color 
transformation results determined to be correct is 
selected as a color transformation result L*a*b*4 

25 (L4,a4,b4) with respect to the RGB signal on the basis 
of the above distances Dl to D3, so that a processing 
time required to create a relationship is shortened. 
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When the number of divided regions is large, 
a time required for the process for calculating color 
transformation results in the creation process unit 
143 increases in proportion thereto. In order to 
5 shorten the processing time, one at the largest 
distance from the boundary surface among the correct 
color transformation results is determined to have 
the highest reliability among the plural correct 
color transformation results, and that color 

10 transformation result is determined as a final color 
transformation result L*a*b*4 (L4,a4,b4). 

When there is no color transformation result 
(L*a*b* values) determined to be correct with respect 
to one RGB signal (grid point signal values) , a color 

15 transformation result L*a*b*4 (L4,a4,b4) with 
respect to the RGB signal is calculated on the basis 
of reciprocals of the distances Dl to D2 between the 
color transformation results obtained with respect 
to the RGB signal and the boundaries (boundary 

20 surfaces) to which the respective color 
transformation results belong. Even if all the color 
transformation results are incorrect, it is possible 
to calculate a color transformation result L*a*b*4 
(L4,a4,b4) having higher reliability . Namely, when 

25 all the color transformation results are incorrect, 
one at the shortest distance from its boundary surface 
Bl, B2 or B3 among the incorrect color transformation 
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results is determined to have higher reliability. 

When there is no color transformation result 
(L*a*b* values) determined to be correct with respect 
to one RGB signal (grid-like signal values) , one among 
5 color transformation results is selected as a color 
transformation result L*a*b*4 (L4,a4,b4) with 
respect to the RGB signal on the basis of reciprocals 
of the distances Dl to D3 . Even if all the color 
transformation results are incorrect, it is possible 
10 to obtain a color transformation result L*a*b*4 
(L4,a4,b4) having higher reliability within a short 
time. In such case, a color transformation result 
closest to the boundary surface is determined to have 
the highest reliability, and a color transformation 
15 result having the largest reciprocal of the distance 
is selected, as stated above. 
[4] Description of modifications of the embodiment 
Next, description will be made of 
modifications of the embodiment with reference to 
20 Figs. 13 through 17. 

[4-1] Description of a first modification 

In the above embodiment, three regions 1 to 
3 obtained by dividing the whole L*a*b* space 
. according to hue angle are beforehand defined. A 
25 first modification of this embodiment will be 
described by way of example where the whole L*a*b* 
space is divided into two regions 1 and 2 according 
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to chroma, with reference to Figs. 13 through 16. 
Fig. 13 is a diagram for illustrating an example of 
region division according to chroma. Figs. 14 
through 16 are diagrams for illustrating practical 
5 examples of an interpolating process according to 
the first modification of the embodiment. 

According to the first modification, two 
regions 1 and 2 obtained by dividing the whole L*a*b* 
space are beforehand defined and set (dividing step) . 

10 As shown in, for example, Fig. 13, the L*a*b* space 
is divided by boundaries (practically, boundary 
surfaces in parallel to the L* axis) Bll and B12 
according to chroma C* [ = ( a * 2 + b* 2 ) 1/2 ] to obtain the 
two regions 1 and 2 . 

15 As shown in Fig. 13, the region 1 is a region 

whose chroma C* is not less than a predetermined value 
C*l (for example, 10.0) in the L*a*b* space. The 
region 2 is a region whose chroma C* is less than 
a predetermined value C*2 (for example, 20.0) in the 

20 L*a*b* space. Therefore, a region having a width 
C*2-C*l (10.0) overlapping on each other is added 
to and formed in the regions 1 and 2 shown in Fig. 
13 . 

Even when the regions 1 and 2 obtained by 
25 dividing the whole L*a*b* space according to chroma 
C* is used, the color transformation table 22 for 
the image scanner 60 is created basically in the 
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similar procedures (steps Sll to S25 in Fig. 5) to 
the above by the color transformation table creating 
apparatus 1 described above with reference to Figs. 

I through 12. 

5 Hereinafter, procedures for creating the 

color transformation table 22 will be described 
briefly . 

[4-1-1] Description of procedures for spectral 
characteristics estimation 
10 The spectral characteristics estimation unit 

II (step S16) estimates the spectral characteristics 
on the basis of several tens of color signals and 
several tens of kinds of spectral reflectance in the 
similar procedures to the above in [2-2] [refer to 

15 the Expression (10)]. 

[4-1-2] Description of procedures for spectral 
reflectance classification 

The spectral reflectance classification unit 
12 (step S18) in the first modification which divides 

20 the space into regions as shown in Fig. 13 determines 
which region 1 or 2 shown in Fig. 13 L*a*b* values 
of each color patch belong to, and classifies spectral 
reflectance of each color patch into either one of 
two groups 1 and 2 according to its region. 

25 Namely, spectral reflectance whose L*a*b* 

values are determined to belong to the region 1 (C* 
^C*l) is classified into the group 1. Spectral 
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reflectance whose L*a*b* values belong to the region 
2 (C*<C*2) is classified into the group 2. 

[4-1-3] Description of procedures for color 
transformation formula creation and procedures for 
relationship creation 

The relationship creation unit 141 (step S21) 
calculates L*a*b* values corresponding to grid-like 
RGB values in the almost similar procedures to the 
above in [2-4] . 

In the first modification, the relationship 
1 of the whole space is calculated on the basis of 
grid-like RGB values, spectral characteristics of 
the image scanner 60 and principal component vectors 
calculated from spectral reflectance of the group 
1. Similarly, the relationship 2 of the whole space 
is calculated on the basis of grid-like RGB values, 
spectral characteristics of the image scanner 60 and 
principal component vectors calculated from spectral 
reflectance of the group 2. 

Whereby, two sets of L*a*b* values (L*a*b:l, 
L*a*b*2) are calculated as a result of color 
transformation of one set of grid-like RGB values 
using two color transformation formulas. 

Procedures for calculating a relationship of 
one group (spectral reflectance group) is as 
described above with reference to the flowchart 
(steps S41 to S47) shown in Fig. 7, detailed 
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descriptions of which are thus omitted. 

[4-1-4] Description of procedures for 
determination and procedures for a creating process 
The creating process unit 143 (step S24) , for 
5 example, selects either one having higher reliability 
between two sets of L*a*b* values created with respect 
to one set of RGB values by the relationship creation 
unit 141 to determine one set of L*a*b* values. 

Prior to the selecting process, the 
10 determination unit 142 (step S23) determines whether 
a result of color transformation (L*a*b* values) of 
the RGB signal obtained by the relationship creation 
unit 141 belongs to a region corresponding to its 
color transformation formula, in the first 
15 modification. When the result belongs to the 
corresponding region, it is determined that the 
result of color transformation is correct. 

Assuming here that a result of color 
transformation of the relationship 1 is L*a*b*l 
20 (Ll,al,bl), a distance between L*a*b*l (Ll,al,bl) 
and the boundary surface Bll of the region 1 is Dl, 
and, similarly, a result of color transformation of 
the relationship 2 is L*a*b*2 (L2,a2,b2) and a 
distance between L*a*b*2 (L2,a2,b2) and the boundary 
25 surface B12 of the region 2 is D2 , as shown in Figs. 
14 through 16, for example. 

The creation process unit 143 selects and 
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determines one result of color transformation on the 
basis of a result of determination by the 
determination unit 142 and the distances Dl and D2, 
in the following manner. 
5 When two results of color transformation 

L*a*b*l (Ll,al,bl) andL*a*b*2 (L2,a2,b2) belong to 
the regions 1 and 2, respectively, and the two results 
are determined to be correct, it is determined that 
one at a larger distance between the distances Dl 

10 and D2 from the respective boundaries Bll and B12 
has higher reliability. Namely, the process is 
performed according to two cases below. 

® In the case where D1>D2: Color 
transformation result L*a*b*l (Ll,al,bl) is 

15 selected. 

(2) In the case where DKD2: Color 
transformation result L*a*b*2 (L2,a2,b2) is 
selected . 

The example shown in Fig. 14 corresponds to 
20 the above case (2) since D2>D1 although two results 
of color transformation L*a*b*l (Ll,al,bl) and 
L*a*b*2 (L2,a2,b2) belong to the corresponding 
regions 1 and 2, respectively, so that the color 
transformation result L*a*b*2 (L2,a2,b2) is 
25 selected. 

When D1-D2, an interpolation equation as 
shown by Expression (18) above, for example, is 
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employed, and an average value of the two color 
transformation results L*a*b*l (Ll,al,bl) and 
L*a*b*2 (L2,a2,b2) is calculated . The average value 
is then selected as a final color transformation 
5 result. 

When either one of the color transformation 
results L*a*b*l (Ll,al,bl) and L*a*b*2 (L2,a2,b2) 
is determined to be correct, the color transformation 
result determined to be correct is selected as it 
10 is as a final color transformation result. 

In the example shown in Fig. 15, the color 
transformation result L*a*b*l (Ll,al,bl) is selected 
since only the color transformation result L*a*b*l 
(Ll,al,bl) belongs to the corresponding region 1. 
15 When both of the two color transformation 

results L*a*b*l (Ll,al,bl) and L*a*b*2 (L2,a2,b2) 
do not belong to the corresponding regions 1 and 2, 
and both of the two are determined to be incorrect, 
it is determined that one at a smaller distance between 
20 the distances Dl and D2 from respective boundaries 
Bll, B12 has higher reliability . Namely, theprocess 
is performed according to two cases below. 

(D In the case where DKD2 (that is, when 
1/D1>1/D2): color transformation result L*a*b*l 
25 (Ll,al,bl) is selected. 

<D In the case where Dl>D2 (that is, when 
1/DK1/D2): color transformation result L*a*b*2 
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(L2,a2,b2) is selected. 

The example shown in Fig. 16 corresponds to 
the above case CD since DKD2 although both of the 
two color transformation results L*a*b*l (Ll,al,bl) 
5 andL*a*b*2 (L2,a2,b2) do not belong to the respect ive 
corresponding regions 1 and 2, thus the color 
transformation result L*a*b*l (Ll,al,bl) is 
selected . 

When D1=D2 (1/D1=1/D2), an average value of 
10 the two color transformation results L*a*b*l 
(Ll,al,bl) andL*a*b*2 (L2,a2,b2) is calculated, for 
example, and the average value is selected as it is 
as a final color transformation result. 

L*a*b* values selected as above are outputted 
15 as a final color transformation result with respect 
to a certain set of RGB values to the memory 20, and 
a relationship with respect to the RGB values is 
determined . 

The above process is performed on all 
20 grid-like RGB values (grid point signal values), 
whereby the color transformation table 22 for the 
image scanner 60 is created in the memory 20. 

[4-1-5] Description of functions and effects of 
the first modification 
25 The above first modification provides 

functions and effects similar to those provided by 
the embodiment described above. Additionally, use 
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of the regions 1 and 2 obtained by dividing the whole 
L*a*b* space according to chroma makes accurate 
calculation of relationships possible, according to 
the first modification. Namely, when principal 

5 component vectors are calculated on the basis of a 
spectral reflectance group obtained by classifying 
spectral reflectance of a color chart into chroma 
regions, a rate of contribution of the principal 
component vectors to the spectral reflectance group 

10 belonging to the corresponding chroma region 1 or 
2 is high, as a matter of course. Therefore, such 
principal component vectors can express a spectral 
reflectance group belonging to the corresponding 
chorma region 1 or 2 more accurately, so that 

15 relationships can be calculated more accurately, as 
a result. Such effect is remarkable when a color 
patch including colors having various kinds of 
chroma . 

[4-2] Description of a second modification 
20 In the embodiment and the first modification 

described above, the regions 1 to 3 obtained by 
dividing the whole L*a*b* space according to hue angle 
or chroma are beforehand defined. In a second 
modification, description will be made by way of an 
25 example where the whole L*a*b* space is divided into 
two regions 1 and 2 according to lightness, with 
reference to Fig. 17. Fig. 17 is a diagram for 
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illustrating an example of region division according 
to lightness in the second modification of this 
embodiment . 

In the second modification, two regions 1 and 
5 2 obtained by dividing the whole L*a*b* space are 
beforehand defined and set, similarly to the first 
modification (dividing step) . These two regions 1 
and 2 are obtained by dividing the L*a*b* space by 
boundaries (practically boundary surfaces in 

10 parallel to the a*b* plane) B21 and B22 according 
to lightness, as shown in Fig. 17, for example. 

As shown in Fig. 17, the region 1 is a region 
whose lightness L* is not less than a predetermined 
valueL*lintheL*a*b* space. The region2 isaregion 

15 whose lightness L* is less than a predetermined value 
L*2 in the L*a*b* space. Therefore, a region having 
a width of L*2-L*l overlapped on each other is added 
to and formed in the regions 1 and 2 shown in Fig. 
17 . 

20 When the regions 1 and 2 obtained by dividing 

the whole L*a*b* space according to lightness are 
used, the color transformation table creating 
apparatus 1 described above with reference to Figs. 
1 through 12 creates the color transformation table 

25 22 for the image scanner 60 according to basically 
the similar procedures to those described above 
(steps Sll to S25 in Fig. 5) . That the number of the 
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regions is two is the same as the first modification 
described above. In the second modification, the 
color transformation table 22 is created in the almost 
similar procedures to those in the embodiment and 
5 first modification described above. 

The second modification provides similar 
functions and effects to the above-described 
embodiment. Additionally, use of the regions 1 and 
2 obtained by dividing the whole L*a*b* space 

10 according to lightness makes accurate calculation 
of relationships possible according to the second 
modification. Namely, when principal component 
vectors are calculated on the basis of a spectral 
reflectance group obtained by classifying a color 

15 chart into lightness regions, a rate of contribution 
of the principal component vectors to a spectral 
reflectance group belong to the corresponding 
lightness region 1 or 2 is high, as a matter of course. 
Therefore, such the principal component vectors can 

20 express a spectral reflectance group belonging to 
a corresponding lightness region 1 or 2 more 
accurately, so that a relationship is accurately 
calculated, as a result. Such effect is remarkable 
when a color chart including colors having various 

25 kinds of lightness is used. 

[4-3] Description of a third modification 

In the above embodiment, the color 
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transformation table 22 is created as one correlating 
an RGB signal with L*a*b* values. In a third 
modification, a color transformation table is created, 
which correlates an RGB signal with spectral 
5 reflectance (spectral reflectance corresponding to 
L*a*b* values as a color transformation result) 
calculated using the above Expression (15). 

When a color transformation table correlating 
an RGB signal with spectral reflectance is created, 

10 the color transformation table creating apparatus 
1 described above with reference to Figs. 1 through 
12 creates the color transformation table for the 
image scanner 60 according to basically the similar 
procedures to the above (steps Sll to S25 in Fig. 

15 5) . 

Note that color transformation formulas used 
in the relationship creation unit 141 (step S21) are 
the above Expressions (14) and (15) . Since the 
determination unit 142 (step S23) requires L*a*b* 

20 values corresponding to RGB values in order to 
determine correctness of spectral reflectance, the 
relationship creation unit 141 or the determination 
unit 142 calculates L*a*b* values corresponding to 
RGB values as a color transformation result in the 

25 third modification, similarly to the above-described 
embodiment . 

When the creation process unit 143 (step S24) 
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performs the interpolating operation process 
(interpolating process) , Expressions (21) to (24) 
below are used instead of the above Expressions (17) 
to (20) . 

5 Assuming here that three kinds of spectral 

reflectance RE Fl (REF1 380 , REFI390/ REFI730), REF2 
(REF2 38 0 / REF2 390 , REF2 730 ) / andREF3 (REF3 38 o, REF3 390 , 

REF3 730 ) are obtained as three relationships 1 to 
3 with respect to one set of RGB values. L*a*b* values 
10 (color transformation results) corresponding to 
these kinds of spectral reflectance are assumed to 
beL*a*b*l (Ll,al,bl) , L*a*b*2 (L2,a2,b2) andL*a*b*3 
(L3 , a3 , b3 ) . 

Similarly to the examples shown in Figs. 9 
15 through 12, a distance between L*a*b*l (Ll,al,bl) 
and the boundary surface Bl of the region 1 is Dl, 
a distance between L*a*b*2 (L2,a2,b2) and the 
boundary surface B2 of the region 2 is D2, and a 
distance between L*a*b^3 (L3,a3,b3) and the boundary 
20 surface B3 of the region 3 is D3. 

When L*a*bn <Ll,al,bl), L*a*b*2 (L2, a2,b2) 
and L*a*b*3 (L3,a3,b3) belong to the respective 
regions 1 to 3, as shown in, for example, Fig. 9, 
and L*a*b*l (Ll,al f bl), L*a*b*2 (L2,a2,b2) and 
25 L*a*b*3 (L3,a3,b3) are determined to be correct , one 
kind of spectral reflectance REF4 (REF4 3 80 r REF4 3 go, 
REF4730) is calculated and determined on the basis 
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of the distances Dl to D3 by using an interpolation 
formula as shown by Expression (21) below. 

REF4 j— ( REF1 j * D14-REF2 j * D2 + REF3 j *D3 ) / (D1 + D2 + 
D3) ... (21) 
5 (where j = 380, 390, 730.) 

When L*a*b*l (Ll,al,bl) and L*a*b*2 
(L2,a2,b2) belong to the corresponding regions 1 and 
2, respectively, as shown in, for example, Fig. 10, 
and L*a*b*l (Ll,al,bl) and L*a*b*2 (L2,a2,b2) are 
10 determined to be correct, one kind of spectral 
reflectance RE F4 ( RE F4 38 q , REF4 390 , REF4 730 ) is 

calculated and determined on the basis of the 
distances Dl and D2 by using an interpolation formula 
as shown by Expression (22) below. 
15 REF4 j= (REFl j *Dl+REF2 j *D2) / (D1 + D2) ... (22) 

(where j = 380, 390, 730.) 

When only L*a*b*l (Ll,al,bl) belongs to the 
corresponding region 1 as shown in, for example, Fig. 
11, and only L*a*b*l (Ll,al,bl) is determined to be 
20 correct, REF1 ( REFl 38 o / REF1 390 , REFI730) is 

determined as it is as REF4 (REF4 3 80rREF4 3 9o, REF4 730 ) 
as shown by Expression (23) below. 

REF4 j=REFl j ... (23) 
(where j - 380, 390, 730) 
25 When none of L*a*b*l (Ll,al,bl), L*a*b*2 

(L2,a2,b2) and L*a*b^3 (L3,a3,b3) belong to the 
corresponding regions 1 to 3, respectively, and all 
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of L*a*b*l (Ll,al,bl), L*a*b*2 (L2 / a2 / a3) andL*a*b*3 
(L3,a3,b3) are determined to be incorrect, one kind 
of spectral reflectance REF4 (REF4 38 o, REF4 390 , 
REF4730) is calculated and determined on the basis 
5 of reciprocals of the distances Dl to D3 by using 
an interpolation formula as shown Expression (24) 
below . 

REF4j= (REFlj* (1/D1) +REF2j* (l/D2)+REF3j*(l/ 
D3) ) / ( (1/D1) + (1/D2) + (1/D3) ) ... (24) 
10 (where j = 380, 390 , 730.) 

Spectral reflectance REF4 (REF4 380 , REF4 390 , 
REF4 730 ) obtained in the above interpolating process 
is outputted as final spectral reflectance with 
respect to one set of RGB values to the memory 20, 
15 and a relationship with respect to the RGB values 
is determined. 

The above process is performed on all 
grid-like RGB values (grid point signal values), 
whereby a color transformation table correlating RGB 
20 values with spectral reflectance is created in the 
memory 2 0 . 

When there are a plurality of color 
transformation results (L*a*b* values) determined 
to be correct, one kind of spectral reflectance REF4 
25 (REF4 380 , REF4 390 , REF4 7 3o) is calculated and 

determined on the basis of the distances Dl to D3 
as shown by Expression (21) or (22), here. 
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Alternatively, it is possible to determine that a 
color transformation result at the largest distance 
(color transformation result farthest away from the 
boundary) has the highest reliability, and select 
and determine it as the spectral reflectance REF4 
(REF4 38 o, REF4 390 , REF4 73 o). 

In the example shown in Fig. 9, REF1 (REF1 380 , 
REF1 390 , REFI730) is selected and determined as final 
spectral reflectance REF4 (REF4 380 , REF4 390 , 
REF4 730 ) since D1>D3>D2 in the above selecting process . 
In the example shown in Fig. 10, REF1 (REF1 3S0 , REF1 390 , 
REFI730) is selected and determined as final 
spectral reflectance REF4 (REF4 380 , REF4 390 , 
REF4 73 o) since D1>D2. 

Similarly, when there is no color 
transformation result determined to be correct, one 
kind of spectral reflectance REF4 (REF4 380 , REF4 390 ; 

REF4 7 3o) is calculated and determined on the basis 
of reciprocals of the distances Dl to D3 as shown 
by Expression (20) . Alternatively, it is possible 
to determine a color transformation result having 
the largest reciprocal of the distance (that is, a 
color transformation result at the shortest distance 
from and closest to the boundary) to be one having 
the highest reliability, and select and determine 
it as the spectral reflectance REF4 (REF4 380 , REF4 390 , 
REF4 730 ) . 
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In the example shown in Fig . 12, REF3 (REF3 380 / 
REF3390/ REF3 730 ) is selected and determined as final 
spectral reflectance REF4 (REF4 3a0 r REF4 390 , 
REF4 730 ) in the above selecting process since 
5 D3<DKD2. 

The third modification provides the similar 
functions and effects to the above-described 
embodiment, as above. Additionally , according to 
the third modification, a color transformation table 

10 correlating RGB signals with spectral reflectance 
REF4 (REF4 380 , REF4 390 , REF4 7 3o) as a color 

transformation result is created, so that the color 
transformation result is obtained as spectral 
reflectance. It is therefore possible to substitute 

15 the image scanner 60 for a simplified colorimeter. 

When there are a plurality of color 
transformation results (L*a*b* values) determined 
to be correct with respect to one set of RGB values 
(grid point signal values) , spectral reflectance REF4 

20 (REF4 38 o/ REF4 390 , REF4 730 ) of the RGB signal is 

calculated on the basis of the distances Dl to D3 
between the two or three color transformation results 
determined to be correct and boundaries (boundary 
surfaces) of regions to which the respective color 

25 transformation results belong, whereby gradation of 
the spectral reflectance is assured. 

When there are a plurality of color 



94 



transformation results (L*a*b* values) with respect 
to one set of RGB values (grid point signal values) 
determined to be correct, one of plural kinds of 
spectral reflectance according to the plural color 
5 transformation results determined to be correct is 
selected as the spectral reflectance REF4 (REF4 380 , 
REF4 390 , REF4 7 3o) on the basis of the distances Dl 
to D3 between the plural color transformation results 
determined to be correct and boundaries (boundary 

10 surfaces) of regions to which the respective color 
transformation results belong, whereby a processing 
time required for relationship creation is shortened . 

When there is no color transformation result 
(L*a*b* values) with respect to one set of RGB signals 

15 (grid point signal values) determined to be correct, 
spectral reflectance REF4 (REF4 380 , REF4 390 , 
REF4 73 o) of the RGB signal is calculated on the basis 
of reciprocals of the distances Dl to D3 between a 
plurality of color transformation results obtained 

20 with respect to the RGB signal and boundaries 
(boundary surfaces) of regions to which the respective 
color transformation results belong. Even if all the 
color transformation results are incorrect, it is 
possible to calculate the spectral reflectance REF4 

25 (REF4 380 , REF4 390 , REF4 730 ) having higher 

reliability . 

When there is no color transformation result 
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(L*a*b* values) with respect to one RGB signal (grid 
point signal values) determined to be correct, one 
of plural kinds of spectral reflectance corresponding 
to the plural color transformation results is 
5 selected as the spectral reflectance REF4 (REF4 380 , 
REF4 390 , REF4 730 ) with respect to the color signal 

on the basis of reciprocals of the distances Dl to 
D3 between the plural color transformation result 
obtained with respect to the color signal and 

10 boundaries (boundary surfaces) of regions to which 
the respective color transformation results belong. 
Even if all the color transformation results are 
incorrect, it is possible to obtain the spectral 
reflectance REF4 (REF4 38 o, REF4 390 , REF4 730 ) having 

15 higher reliability within a short time. 

[4-4] Description of a fourth modification 

In the embodiment and the first to third 
modifications described above, color transformation 
results (L*a*b* values ) and distances from boundaries 

20 of regions are used when the creation process unit 
143 performs the interpolating process or the 
selecting pro cess on the color transformation results 
or spectral reflectance. Instead of the distance, 
a value of A E94 obtained by transforming the distance 

25 may be used. 

In such case, a color transformation table 
that attains more accurate color transformation when 



96 



viewed from man can be created. AE04 has a closer 
color difference when viewed by man than a distance 
( A Eab) in the L*a*b* space. For this, a color 
transformation result using AE94 is more accurate 
5 when viewed by man. 
[4-5] Others 

The present invention is not limited to the 
above examples, but may be modified in various ways 
without departing from the scope of the invention. 
10 In the above-described embodiment, the 

transformation source color space is the RGB space 
and the transformation target color space is the 
L*a*b* space, for example. However, the present 
invention is not limited to this. The present 

15 invention may be applied when a color transformation 
table correlating color signals in various color 
spaces with each other is created. In such case, the 
present invention can provide the similar functions 
and effects to the above-described embodiment. 

20 In the above embodiment, the transformation 

target color space is divided into three or two. 
However, the present invention is not limited to this . 
It is possible to divide the color space into four 
or more regions. The above embodiment has been 

25 described by way of example where the transformation 
target color space is divided according to only hue 
angle, only chroma and only lightness. However, it 
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is alternatively possible to combine at least two 
parameters among hue angle, chroma and lightness, 
and divide the transformation target color space into 
a plurality of regions according to these parameters. 
5 The above embodiment has been described by 

way of example where principal component vectors, 
spectral characteristics or spectral reflectance are 
36 discrete data obtained correspondingly to 
wavelengths (380, 390, 730 nm) for every 10 nm 

10 in the visible radiation wavelength range 380-730 
nm. The present invention does not limit principal 
component vectors, spectral characteristics and 
spectral reflectance to discrete data, as above. 

In the above embodiment, the determination 

15 unit 142 determines using L*a*b* values calculated 
by the relationship creation unit 141. 
Alternatively, the determination unit 142 may 
calculate L*a*b* values. 
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